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I  INTRODUCTION 


The  goal  of  this  project  is  the  development  of  optical  signal  processing  architectures 
and  techniques  that  are  suitable  for  processing  information  in  two  dimensions  and  yet 
utilize  principally  one  dimensional  transducers.  The  method  that  we  use  to  accomplish  this 
is  time  and  space  integrating  processing  and  it  allows  the  implementation  of  very  powerful 
analog  signal  processing  architectures  using  well  developed  acoustooptic  devices  as  input 
electronic-to-optical  transducers.  This  in  turn  makes  it  possible  in  many  instances  to  apply 
in  the  immediate  future  these  architectures  to  applications  that  require  two  dimensional 
signal  processing  capability.  We  have  applied  our  methods  to  synthetic  aperture  radar, 
image  recognition,  two  dimensional  spectrum  analysis  of  one  dimensional  signals  and 
adaptive  phased  arrays.  In  previous  years  we  have  developed  several  specific  acoustooptic 
architectures  for  each  application  area  and  we  have  experimentally  demonstrated  at  least 
one  architecture  in  each  category.  A  direct  consequence  of  our  work  in  this  area  has 
been  the  initiation  of  development  programs  in  industrial  and  government  laboratories  to 
build  architectures  that  we  have  developed  under  this  contract  for  specific  applications. 
For  instance,  there  are  at  least  five  acoustooptic  synthetic  aperture  radar  processors  (a 
system  developed  at  Caltech  under  this  contract)  that  have  been  constructed  in  industrial 
laboratories  and  there  are  active  development  programs  in  each  of  the  four  application  areas 
we  mentioned  above.  In  this  report  we  present  our  results  in  each  of  the  four  application 
areas  during  the  past  year. 


The  emphasis  of  our  work  in  the  radar  imaging  area  has  shifted  to  inverse  synthetic 
aperture  radar  (ISAR).  In  ISAR  the  object  that  is  being  imaged  is  moving  usually  in 
a  trajectory  that  is  not  completely  known.  Therefore  ISAR  is  a  much  more  demanding 
signal  processing  problem  than  SAR  because  the  trajectory  of  the  moving  object  needs  to 
be  estimated  in  addition  to  forming  the  radar  image  once  the  trajectory  is  known.  The 
approach  that  we  have  developed  is  to  cast  the  ISAR  problem  as  an  optimization  problem. 
An  image  is  formed  using  the  current  estimate  of  the  trajectory,  and  then  an  “energy’’ 
function  is  calculated  that  is  mininum  when  the  image  is  sharply  focused.  The  estimate  of 
the  trajectory  is  then  modified  to  decrease  the  energy  function  and  thus  focus  the  image. 
Depending  on  the  characteristics  of  the  target  trajectory,  it  can  take  hundreds  of  iterations 
for  the  algorithm  to  converge.  In  SAR,  where  the  trajectory  is  assumed  known,  only  a 
single  iteration  is  required.  Therefore  the  computational  speed  that  is  required  for  ISAR 
(at  least  when  our  algorithm  is  used)  is  several  orders  of  magnitude  higher  than  SAR. 
Current  digital  signal  processing  techniques  have  sufficient  processing  power  to  do  the 
SAR  problem  in  real  time.  The  reasons  that  optics  can  be  advantageous  in  certain  SAR 
applications  is  not  speed,  but  rather  low  power  consumption  and  size.  For  real  time  ISAR 
imaging  however  in  applications  where  the  target  trajectory  is  not  known  the  required 
speed  is  beyond  the  capabilities  of  digital  signal  processing  techniques.  An  analog  optical 
implementation  however  has  sufficient  processing  power  to  tackle  the  problem  in  real  time 
with  reasonably  low  power  requirement  and  system  complexity.  The  energy  formulation 
of  the  ISAR  problem  is  described  in  section  If. 


During  the  past  year  we  completed  the  experimental  demonstration  of  the 
programmable  acoustooptic  image  correlator.  In  this  system  an  input  image  is  detected  by 
a  TV  camera  and  applied  to  the  optical  system  as  a  video  signal  through  an  acoustooptic 
device.  The  reference  is  stored  in  a  digital  memory  and  entered  into  the  system  through  a 
one  dimensional  array  of  LEDs.  In  section  III  we  show  that  if  the  reference  is  stored  with 
only  one  bit  of  accuracy,  then  the  performance  degrades  only  marginally.  The  capacity  of 
the  correlator  (defined  as  the  maximum  number  of  distinct  objects  it  can  classify  correctly) 
is  reduced  by  a  factor  of  only  .on.  In  addition  the  traditional  training  algorithms  that  exist 
for  programming  the  correlator  can  be  readily  modified  to  accommodate  the  binary  format 
of  the  reference.  In  section  III  we  present  a  binary  version  of  the  perceptron  algorithm.  The 
significance  of  the  above  results  is  that  the  design  of  the  electronic  portion  of  the  optical 
correlator  is  simplified  drastically.  The  experimental  demonsrtration  of  the  acoustooptic 
correlator  performing  multiple  target  recognition  is  also  reported  in  section  III. 

Two  different  acoustooptic  architectures  for  performing  two  dimensional  spectrum 
analysis  were  demonstrated.  The  first  is  a  Mach-Zender  interferometer  with  two 
orthogonally  oriented  acoustooptic  devices  in  each  of  its  arms.  The  two  arms  of  the 
interferometer  interfere  on  a  two  dimensional  CCD  where  the  2-D  Fourier  transform 
accumulates  through  temporal  integration.  The  description  of  this  system  and  its 
experimental  demonstration  has  been  published  as  a  paper  (reference  2  in  section  IV). 
The  second  architecture  is  a  multiplicative  one  (the  light  from  one  acoustooptic  device  is 
then  rediffracted  from  a  second  that  is  orthogonally  oriented.  This  second  architecture 
is  described  in  section  IV  and  it  is  more  insensitive  to  mechanical  vibrations  than  the 
interferometric  approach.  The  interferometer  however  has  its  own  advantageous  features 
which  are  better  light  efficiency  and  simpler  demodulation  requirements  at  the  output  to 
separate  the  signal  from  teh  bias  terms. 

Section  V  reports  our  progress  in  the  area  of  adaptive  filters.  In  previous  years  we  have 
demonstrated  single  channel  acoustooptic  filters  that  utilized  a  combination  of  temporal 
and  spatial  integration.  This  year  we  concentrated  on  broadband  phased  arrays.  In  this 
case  jammers  may  be  present  anywhere  within  a  bandwidth  as  high  as  one  gigahertz. 
The  adaptive  processor  must  then  place  nulls  in  both  the  temporal  frequency  domain  and 
the  direction  of  arrival  of  the  jammer  to  effectively  cancel  the  interference.  Again  this 
is  a  problem  that  is  extremely  demanding  computationally  and  analog  optics  provides  a 
uniquely  effective  solution.  For  example,  for  a  GHz  bandwidth.  100  antenna  elements  in 
the  array,  and  100  taps  within  each  channel,  the  computational  rate  required  simply  to  do 
the  broadband  beamforming  is  101’’  multiplications  per  second.  In  addition  the  weights 
for  each  weight  must  be  adaptively  calculated  in  real  time  which  further  increases  the 
computational  load.  It  is  very  unlikely  that  a  technology  other  than  analog  acoustooptic 
signal  processors  will  be  capable  in  the  forseeable  future  of  achieving  this  processing  rate 
with  reasonable  power  consumption  and  system  complexity. 


II  ISAR  IMAGING 


II. 0  Introduction 

The  theory  and  technology  of  producing  synthetic  aperture  radar  images  of  stationary 
targets  from  moving  platforms  have  been  well  developed  and  utilized1.  As  SAR  imaging 
requires  precise  knowledge  of  the  relative  motion  between  radar  and  target  it  is  not  possible 
to  directly  produce  SAR  images  of  targets  which  themselves  move  beyond  the  control  of 
the  imaging  system.  Inverse  synthetic  aperture  radar  (ISAR)  is  an  attempt  to  overcome 
this  difficulty2,3. 

II.  1  SAR  Imaging  Fundamentals 


Consider  a  three-dimensional  rigid  target  whose  reflectivity  at  any  point  in  space  is  given, 
in  a  coordinate  system  fixed  on  the  target,  by  /(r)  (That  /(r)  does  not  depend  on  time  is 
a  statment  of  the  constraint  that  the  target  is  a  rigid  body).  Any  relative  motion  between 
this  target  and  an  imaging  radar  will  be  described  as  seen  in  the  target’s  coordinate  system; 
that  is  the  radar  moves  along  a  trajectory,  R(t),  while  the  target  remains  fixed.  It  will  be 
assumed  that  the  change  in  R(t)  over  the  duration  of  a  single  pulse  can  be  neglected 

If  the  radar  transmits  a  pure  tone,  g[t )  =  ej2’rt'"t,  and  the  target  consists  of  a  small 
volume  d3r  centered  on  r  with  reflectivity  f(r)  the  received  signal  will  be  given  by 

-s(t)  -  d3r/(rje; 27r""(f  2|r' 

If  the  target  consists  of  a  collection  of  such  volumes,  and  if  multiple  reflections  and  shad¬ 
owing  effects  are  neglected,  then  the  received  signal  will  be  given  by 

s (t)  -  J  d*r  f(r)ej2iru"(t  2  r'  *!/,:). 

If  R  3>  r  then  the  far- field  approximation  gives  |r*—  R\  s;  R  -  eR  •  r  where  e‘n  is  the  unit 
vector  in  the  R  direction.  Consequently  the  received  signal  is  given  by 


In  other  words  the  radar  receives  the  (spatial)  Fourier  transform  of  the  object,  at  spatial 
frequency  2u"  2  along  the  ray  in  Fourier  space  defined  by  e'/?.  on  a  temporal  carrier. 


If  the  radar  instead  transmits  a  signal  with  a  finite  bandwidth, 

■  H  ’2 


/n 

H,2 


<lv  G(i/)e 


, )  t 


it  follows  that  the  received  signal  will  be  given  by 


s(t)  =  j  du'c{u)e-}2wu2R/cF('^-^-l^hR)\e}2K‘/t. 

J  B/2  c 

Mixing  s(t)  with  e-;2,ri/"(f~2fl/‘;).  Fourier  transforming  and  multipling  the  result  by 

*(„)  =  /  “»/*  "  <  B/2 ; 

1 0,  otherwise 

will  produce 

W(i/)  =  F(— ^^e-ft)rect(^). 

If  the  radar  does  this  for  N  different  rays,  0  <  k  <  N  —  1,  it  will  produce  samples 
of  the  target’s  Fourier  transform  over  some  volume  of  Fourier  space.  These  can  then  be 
inverse  transformed  producing  the  function  /(r)  to  within  the  resolution  permitted  by  the 
finite  size  of  the  volume.  This  is  the  basic  principle  of  synthetic  aperture  radar. 

II. 2  Two-Dimensional  Imaging 


The  special  case  of  a  two-dimensional  geometry  is  considered  in  this  paper  because  —  I) 
It  has  important  practical  applications;  almost  all  work  in  the  literature  deals  with  this 
case,  and  2)  it  greatly  reduces  the  complexity  of  the  problem  for  analysis  and  computer 
simulation,  yet  the  results  readily  generalize  to  three-dimensions.  For  a  two-dimensional 
-  cos  sin  9k  j.  Consequently  a  two-dimensional  target,  /(j,y),  gives  rise 
to  a  sequence  of  N  Fourier  samples  described  by 

F(  -  —  -  °  ;cos  0*;,sin  0fc|)rect(i//#) 
c 


rect (is/B)  J  dx  J  dy  f(x,y)ej2*  ‘  (zc 


+  y  ^ i n  ) 


0  <k<  N  1. 


A  very  important  class  of  problems  results  when  the  small  angle  conditions, 

2u 

cos  Bk  =  1,  sin  0k  as  9k,  (  y6)m.lT  -  <  1 

c 

hold.  Then  the  Fourier  samples  “decouple’’  into  a  cascade  of  one  dimensional  transforms 


J  dye12*  yflk  J  dr  rect(/v/ H)e’ 2*  zf(x,y),  0  <  k 


<  N  1. 


Resolution  in  x  is  achieved  via  “range  compression'1 


dve~j2* 


2  ( is  > 


/(x,y)  *  B  sinc( 


x 

Vj2B 


\ 

)■ 


What  is  left  are  the  samples 

=  J  dy  e32*^  ydk  f{x,y)  *  Bs\nc{—^). 

The  quantity  c/25  =  Ax,  the  x  resolution,  defines  the  width  of  a  “range  bin”.  The  rest  of 
the  discussion  will  center  on  achieving  y  resolution  within  any  given  range  bin  (fixed  x). 
For  clarity,  therefore,  explicit  dependence  on  x  will  no  longer  be  shown. 

Resolution  in  y  is  thus  achieved  by  inverting 

Fk  =  J  dy  f{y)e]^-y9k,  0  <  k  <  N  -  1, 

which  is  done  by  simply  inverse  transforming  via 


(1 1 .2.1) 


1 

N 
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If  the  problem  is  well  posed  (if  it  is  at  ail  possible  to  form  an  image)  then  the  quantity 
in  brackets  will  be  a  sinc-like  function  with  some  main  lobe  width  Ay.  Since  the  actual 
function  depends  on  the  Ok  in  a  very  complicated  way,  the  following  approximation  will  be 
employed 


1 

N 


»’"*  *  %  -  y). 


[11.2.2] 
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It  should  always  be  remembered  that  this  is  a  “delta”  function  only  on  a  scale  larger  than 
Ay  (it  really  has  a  finite  width).  In  view  of  (II. 2. 2)  it  follows  that  to  within  the  finite 
resolution  Ay.  /(y)  =s  /(y). 


Notice  the  implicit  assumption  that  R  for  each  pulse  —  Rk  ---  (needed  for  mixing 
and  filtering)  and  the  $k  are  known.  This  is  the  case  in  SAR  imaging.  For  ISAR  imaging, 
however,  the  trajectory  is  unknown  a  priori.  What  one  has  instead  are  estimates  /?*.,  §k- 
Since  radar  is  a  range  measuring  device,  accurate  “gross-range”  estimates,  Rk,  are  much 
more  straight  forward  to  obtain  than  are  the  “aspect”  angles,  §k,  through  techniques  such 
as  point-target  referencing  or  echo  correlating.  Therefore  discussion  below  is  confined  to 
estimating  the  aspect  angles. 


5 


The  N  aspect  angle  estimates,  0  —  (Oo, . . . ,  9 s -  1 ),  define  an  N  dimensional  state  space. 
With  each  point  in  state  space  is  associated  an  image  — 

f(y,e)  =  I  dr,  f(r,)  =  I  dr,  f(n)h(y,m'6). 

For  6^9,  h  is  in  general  not  a  “delta”  function,  and  consequently  f(y)  ^  /(y).  In  its  most 
general  form,  then,  the  ISAR  problem  is  to  search  state  space  for  the  point  corresponding 
to  the  “true”  image.  (Notice  that,  because  of  the  small  angle  assumptions,  h  has  a  form 

which  is  invariant  under  y  — ►  y/a,6  — 1 •  a6.  Consequently  all  images  associated  with  a  given 
ray  in  state  space  are  scaled  (in  y  only)  versions  of  a  single  image.  This  fact  is  what  makes 
ISAR  as  it  is  currently  described  in  the  literature2,3  possible.  One  merely  assumes  that 
6k  oc  k\  the  actual  constant  of  proportionality  is  not  needed.  Hence  there  is  no  need  for 
parameter  estimation  and  the  problem  reduces  to  an  inverse  Fourier  transform.  However 
this  approach  only  works  for  trajectories  consistent  with  this  assumption.) 


II. 3  State  Space  Statistics 

Assume  the  range  bin  consists  of  N  point-targets  located  at  positions  yn  =  nAy,  i.e. 


w-i 


f{y)  =  5Z  -  nAy) 


(/7.3.1) 


n=0 


where  /  and  0  are  the  magnitude  and  phase  of  the  reflectance.  Then  the  returns  are  given 

by 

Fk  =  fme34”"e]m*k 


where  4>k  —  ^/■/^k-  If  the  target  is  reconstructed  using  trajectory  estimates  Qk  -  $k  +  f-k 


then  the  image  pixels  are  given  by 

J'/'m  „j{rn<t>k -n4>k)  _ 
m  k 

and  hence 


^  '"~5> 
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In? 


k  l 


N  2 


j{m  n).t>ke~](p  n)*ie  jn(rk  U  ) 


m  p 


If  the  ik  are  treated  as  independent,  identically  distributed,  random  variables  with 
density  p(f.)  then  the  expected  value  of  fn  2  can  be  computed  as 

Fj  f„  2  -  *  Y1  Yl  '*v)  n)0‘e  ]iP  }nUk  ',yr 


The  expected  value  of  the  random  term  is 


=  I  l  f  de  p(e)e  >nt\2  =  ,P(rc)|2,  if  k  jt  l 


if  A:  =  / 


or  more  compactly 


=  |/>(n)|2  +  «H(1  _  |p,(n)|2). 


Substituting  (4)  into  (II. 3. 2)  it  follows  that 


E  |/nH  = 


2,_  I PWl* 


j'l'/'m-t/'p) 


^2J2ej{m-n)<t>ke-]{p-n)4>l  (//.3.4) 


,  U  ~  lp2(n)l2)  ^  e^m~p^k 

m  p  k 


By  (II.2.2), 


Ye3{m~p)*k  *  -V6mp  =>Y52*iim-n)*keiip-n)*,  «  N2SmnSpn. 


Using  this  in  equation  (II. 3. 4)  gives  the  final  result  — 

E[I/„!2M/’|P(»)|2  +  (L-J£WL2)  £ /= 


(//.3.5) 


Equation  (II. 3. 5)  says  two  things  —  errors  in  the  trajectory  estimate:  l)degrade  the 
image  pixels  from  an  intensity  of  f2  to  f2\P(n)\2  (the  first  term),  and  2)give  rise  to  a  noise 
floor  (the  second  term).  As  an  example  let 

p{e)  =  — 

\/2no 


and  (II. 3. 5)  becomes 


P{n)  ■- 


/;«  .  Y.&- 


(//.3.5a) 


The  pixels  farthest  from  the  y  0  axis  (n  large)  have  on  average  both  a  higher  noise 
level  and  a  lower  reconstruction  intensity.  Thus  the  image  deteriorates  from  the  extremities 


.00  OY»V 


inward.  This  is  as  expected  since  a  small  error  in  aspect  angle  causes  a  phase  error  which 
grows  proportional  to  y.  To  make  this  statment  quantitative  a  signal  to  noise  ratio  can  be 
defined  as  the  ratio  of  the  first  term  in  (II. 3. 5)  to  the  second  — 


SNR(n) 


fl  |P(»)| 


Eav  1  -  |P(n)i2 


(I  [.3.6) 


where 


-V  f2 

v  Jm 


N 


For  the  p(c)  used  in  the  above  example  (II. 3. 6)  becomes 


SNR{n)  = 


f2 

J  n 


;(n<T)2  -  1 


(1 1. 3.6a) 


which  clearly  displays  the  pixel  deterioration  relative  to  the  noise  floor  with  increasing  n. 
Figures  II. 1  and  II. 2  show  computer  simulation  verification  of  (II. 3. 5a)  with  N  =  16,  a  = 
0.5,  for  a  range  bin  with  point-targets  at  pixels  0  and  7. 


1 1. 4  Energy  Functions 


ISAR  imaging  can  be  treated  as  an  optimization  problem  with  respect  to  the  state  space 
parameters  0fc.  The  central  task  then  becomes  the  formulation  of  an  “energy”  function 
describing  the  quality  of  the  reconstruction.  Without  some  a  priori  assumptions  about 
the  nature  of  the  image  to  be  reconstructed  it  is  impossible  to  say  whether  a  given  recon¬ 
struction  is  “good”  or  “bad”  since  all  images  are  equally  likely.  Thus  image  reconstruction 
schemes  generally  contain  a  substantial  amount  of  a  priori  assumptions  as  to  the  form 
of  the  original  image.  The  more  specific  the  assumptions  the  more  discerning  an  energy 
function  that  can  be  constructed. 


Our  research  has  so  far  concentrated  on  energy  functions  for  an  important  class  of 
targets  namely  binary-random-phase  targets.  These  targets  have  the  property  that  fn  is 
0  or  I  while  Vu  is  arbitrary.  Such  a  target  model  should  be  well  suited  to  ISAR  imaging 
of  aircraft  in  flight  and  ships  at  sea  since  such  targets  consist  of  a  nearly  uniform,  highly 
reflective  object  (metal)  against  a  poorly  reflective  background  (air  or  water).  It  can  then 
be  stated  a  priori  that  a  good  reconstruction  will  have  fn  2  very  nearly  1  or  0  with  no 
other  values  allowed.  An  energy  function  making  use  of  this  knowledge  is 


.v  1 


E 


.V 


y  Kn  where  F,n 


I !  ','V  !"•  «  A 2 


o: 


I  m 


if  fn  - 


(ii.  u; 
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where  a  and  p  are  constants  to  be  chosen.  Clearly  hi  is  non-negative  and  assumes  its 
global  minimum  value  of  0  when  the  binary  image  is  faithfully  reconstructed.  K  assumes 


r*  -v/V  ,*  .  *  ./  **  *’  -  *  -  -  •  *  *  *  «  *  .  '  '  *  ».  •-*  %' 


the  value  L  when  all  pixel  energies  assume  the  value  a.  The  normalization  factors  in  the 
denominators  insure  that  E  is  continuous  as  \fn\ 2  moves  across  the  threshold. 


Using  (II. 3. 5a),  choosing  a  to  be  Eav,  and  then  substituting  into  (II. 4.1)  gives  En 
(1  -  e-(n<TU)p  and  hence 


for  the  energy  of  the  mean  image.  This  function  is  monotonically  increasing  in  o.  and 
hence  has  no  local  minima.  It  is  also  independent  of  the  distribution  of  binary  point- 
tar vts,  depending  only  on  their  total  number  through  Eav. 


11. 5  Simulation  Results 

We  have  run  computer  simulations  using  the  range  bin  model  given  in  (II. 3.1)  with  N  =  16. 

Starting  with  some  trajectory,  0,  we  added  noise  to  produce  an  initial  estimate,  0O,  then 
computed  the  corresponding  initial  image  and  its  energy  as  given  by  (II. 4.1).  From  there 
an  optimization  algorithm  (having  no  knowledge  of  the  true  trajectory)  was  used  to  search 
state  space  for  progressively  “better”  images.  Figure  II. 3  shows  a  range  bin  consisting  of 
8  point  targets.  Figure  II. 4  is  the  image  corresponding  to  an  initial  trajectory  estimate 
while  the  image  in  figure  II. 5  corresponds  to  a  point  in  state  space  eventually  reached  by 
the  algorithm. 

While  the  above  approach  is  completely  general,  it  requires  a  search  over  a  space 
with  X  degrees  of  freedom.  For  large  X  this  may  not  be  practical.  At  the  expense  of 
generality  we  can  reduce  the  degrees  of  freedom  by  parameterizing  the  trajectory  (e.g. 
0k  -  axk  f  a2k2  +  . .  .)  with  M  parameters  where  \f  <C  N.  Our  goal  is  then  to  obtain 
faithful  estimates  —  a( ,  a2  ■ . .  a\j .  Figures  II. 6  and  II. 7  show  results  for  the  case  Ar 
16,  .V/  —  2,  where  the  target  is  that  in  figure  II. 3;  figure  II. 6  corresponding  to  the  initial 
parameter  estimates  and  figure  II. 7  corresponding  to  the  improved  estimates  obtained  by 
the  optimization  algorithm. 

1 1. 6  Optical  Processing 

nThe  optical  implementation  of  the  matched  filtering  process  described  in  equation 
II. 2  1  is  straight  forward.  All  we  must  do  is  to  place  the  samples  in  the  front  focal  plane  of 
a  Fourier  transforming  lens  such  that  the  transverse  position  of  F k  is  proportional  to  0K . 
That  is  we  must  produce  a  transmittance  in  this  plane  given  by 

T(s)  Fkd(x  aOk ) 

k 

where  a  is  some  constant.  Due  to  the  Fourier  transforming  property  of  the  lens  we  get  in 
its  bar  k  focal  plane  a  field  (x  arid  y  are  the  transverse  coordinates  in  the  front  and  back 


focal  planes  respectively) 


-  /  dxT{ x)e-J2"»Wx"f)  = 


-j2nyek(a/\..f) 


where  \Q  is  the  optical  wavelength  and  /  is  the  focal  length  of  the  lens.  Comparing  this 
expression  to  the  right  hand  side  of  equation  II. 2.1  we  see  that  they  are  identical  except 
for  a  different  scale  factor  in  the  exponent.  As  we  have  already  seen  this  merely  scales  the 
image  produced  accordingly. 

The  transmittance  T(x)  can  be  produced  by  putting  the  samples  on  a  carrier  and 
then  feeding  them  into  an  acoustooptic  device  with  sample  F k  entering  the  device  at  a 
time  proportional  to  Ok  which  is  our  best  current  estimate  of  Ok-  When  all  the  samples 
have  entered  the  device  a  short  laser  diode  pulse  effectively  freezes  the  position  of  the 
samples  and  produces  the  image  f(y). 

If  we  let  the  image  fall  on  a  photodetector  array  we  can  read  out  the  pixel  intensities 
serially,  computing  the  energy  En  in  equation  II. 4.1  as  each  pixel  emerges  and  by  keeping 
a  running  sum  we  will  have  computed  the  total  energy  E  when  the  last  pixel  emerges. 
This  energy  value  is  then  feed  back  to  the  circuitry  which  determines  when  the  samples 
will  be  feed  into  the  acoustooptic  device  at  the  next  iteration. 
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figure  II.l  -  mean  of  100  images  made  with  sigma=0.5 
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figure  II. 2  -  equation  II. 3. 5a  with  sigma=0.5 
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figure  II. 3  -  target  (single  range  bin) 
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figure  II. 4  -  initial  image 
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figure  II. 5  -  "optimized"  image 
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figure  I I. 6 


initial  image  (2  parameter  case) 
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figure  II. 7  -  "optimized  image" 


Ill  LED  BIPOLAR  CORRELATOR 


III.O  Introduction 

Ever  since  Vander  Lugt  demonstrated  optical  matched  filtering  in  1964 l,  the  main  stream 
of  optical  pattern  recognition  has  been  correlation  type  image  recognition.  In  addition  to 
the  ease  of  implementation  by  optics,  matched  filtering  type  correlator  is  also  insensitive 
to  any  shift  in  position  of  the  input  pattern.  Besides  using  a  Fourier  transform  filter, 
matched  filtering  can  also  be  implementated  by  using  a  correlation  plane  filter.  The  success 
of  either  scheme,  however,  is  limited  by  the  availability  of  linear  spatial  light  modulators 
with  sufficient  dynamic  range  and  accuracy.  Subsequently,  we  resort  to  the  use  of  a  bipolar 
spatial  light  modulator.  In  this  report,  we  address  the  issues  of  the  feasibility  of  using  an 
image  plane  bipolar  filter  for  correlation  and  the  performance  of  bipolar  filters  versus  linear 
filters  in  terms  of  output  signal  to  noise  ratio.  We  also  investigate  the  possibility  of  using 
a  single  bipolar  filter  for  multiple  patterns  recognition;  some  algorithms  of  so  doing;  and 
the  storage  capacity  of  the  bipolar  filters  corresponding  to  those  algorithms.  In  the  last 
four  sections,  we  describe  an  optical  system  designed  to  implement  the  linear  perceptron 
classification  procedure  using  a  bipolar  filter. 


III.l  Bipolar  correlators 


We  begin  our  investigation  by  first  analysing  the  performance  of  thresholding  the  input 
image2.  The  performance  that  we  adopt  here  is  SNR ,  the  signal  (correlation  peak)  to  noise 
(additive  noise  and  side-  lobe)  ratio.  Let  f(i,j)  be  the  image  to  be  recognized.  Assume 
f[i,j)  to  be  a  discrete  sequence  of  independent,  identically  distributed  Gaussian  random 
variables  with  zero  mean  and  variance  oj.  The  input  image  to  the  correlator  is  the  sum 
of  f(i,j),  the  image  to  be  recognized,  and  n(i,j),  an  additive  noise  image.  Here,  we  also 
assume  n(i,j)  to  be  an  independent,  discrete  sequence  of  Gaussian  random  variables  with 
zero  mean  and  variance  <72.  The  thresholded  image  is  defined  as 


if  f{ij)  +  n(i,j)  >  0; 
if  otherwise. 


(NIAA) 


The  reference  image  h(i.j)  is  obtained  by  thresholding  f(i,j): 


I,  if  f{ij)  >  0; 

1.  otherwise. 


(IIIA.2) 


The  correlation  between  the  thresholded  in  put  and  reference  images  is  then  given  by 


,Y  ,V 

»(*'♦/)  Yi  •  *',>  *  j') 


[1 1 1  AM) 


The  signal  to  noise  ratio,  defined  as 


SNR'  = 


E2[g{  0,0)1 


i'J' 


^0 


can  be  shown  to  be 


SNR'  =  N2 


r  -tan~l  [S N Rtn) 
l  7 r 


i  2 


where  S N Rin  =  — . 


(///.  1 .4) 


(///.1.5) 


In  order  to  see  the  effect  of  thresholding  the  input  image,  we  compare  SNR'  with  the 
signal  to  noise  ratio  that  is  obtained  if  the  input  image  is  not  thresholded.  It  can  be  shown 
that  in  this  case 


SNR  = 


N2SNR?n 
l  +  SNR*' 


{III.  1.6) 


The  two  SNRs  are  plotted  in  Figure  III.  1.1  as  a  function  of  SNR^n.  As  can  readily  be 
observed  from  the  plot,  the  correlator  with  input  image  thresholded  is  only  marginally 
degraded  for  high  input  SNR. 


III. 2  Bipolar  rotation  invariant  filters 

Mathematically,  any  pattern  f{r,0)  can  be  decomposed  into  an  infinite  sum  of  its  circular 
harmonics3, 

OO 

/M)  =  E  (///.2.1a) 

n  —  -oo 

where 

fn{r)  =  J  f(r,0)ejn*d0.  (///.2.16) 

Thus,  if  the  pattern  is  correlated  with  any  component  of  its  circular  harmonics  /n(r)e  jnri, 
the  magnitude  of  the  output  is  rotation  invariant.  Furthermore,  if  any  component  of 
circular  harmonics  is  recorded  on  a  Fourier  transform  hologram,  the  resulting  filter  is  also 
shift  invariant. 

Computer  generated  bipolar  filters  were  generated  from  the  first  two  orders  of  cir¬ 
cular  harmonics  of  the  letter  A  using  the  algorithm  described  by  eq. (III. 2.1).  Numerical 
correlations  of  the  letters  A,B,C,and  D  and  the  binary  CGHs  were  evaluated.  The  result 
is  presented  in  Figure  111.2. 1.  The  intra-class  recognition  using  the  zeroth  order  circular 
harmonics  is  remarkable.  However,  the  inter-class  discrimination  performance  is  barely 
satisfactory.  In  order  to  obtain  good  discrimination  performance,  cross  correlations  using 
bipolar  filters  generated  from  different  circular  harmonics  have  to  be  evaluated.  A  compro¬ 
mise  between  the  inter-  and  intra-  class  discrimination  performance  will  determine  which 
bipolar  filter  to  be  used. 


III. 3  Generalized  bipolar  filters 

In  this  section,  we  examine  the  possibility  of  designing  linear  bipolar  filters  for  recognizing 
multiple  objects.  One  immediate  extension  of  such  a  design,  if  it  can  be  done,  is  to 
incoporate  information  corresponding  to  multiple  versions  of  a  certain  object  in  a  filter 
to  achieve  invariance  (e.g.  rotation  and/or  scale).  We  investigate  two  schemes  that  may 
accomplish  this  goal. 

The  first  scheme  is  a  derivative  of  the  well  known  2-category  perceptron  invented  in 
the  60’s4.  Basically,  a  2-category  linear  perceptron  is  a  machine  consisting  of  a  weight 
pattern  w  trained  by  a  certain  algorithm  utilizing  the  set  of  all  input  patterns  {fm}  such 
that  the  output 


sgn 


We  confine  our  attention  to  one  of  the  simplest  perceptron  training  algorithm,  which  is 
as  follows:  The  set  of  input  patterns  is  arranged  as  a  repetitive  sequence  of  patterns,  i.e. 
{f 1 ,  f 2 ,  •  •  •  ,fM,f 1  •••}.  The  weight  vector  is  trained  by  the  elements  from  the  sequence 
one  at  each  time, 


wfc  +  1  =  wfc  +  afcffc;  w1  =  some  initial  vector, 


(/ / /.3.2a) 


where 


+  1  if  gk  --  sgn[ w 
=  \  -l  if  gk  =  +1  and  fk  $ 
0  otherwise. 


k  -ffcl  =  -1  and  ik  e 


(III. 3:2b) 


The  algorithm  has  converged  when  the  correct  response  is  obtained  for  each  input  vector. 

A  bipolar  perceptron  can  be  defined  by  modifying  the  linear  perceptron  simply  by 
replacing  w  by  a  bipolar  weight  pattern.  The  training  algorithm  can  also  be  modified  to 


wfc+ 1  -  sgn\ wfc 


krk: 


a*f 


{ill::,  .3) 


The  result  of  the  rate  of  convergence  of  a  bipolar  perceptron  trained  using  the  above 
algorithm  is  shown  in  Figure  III.3.1.  An  ordinate  reading  of  104  iterations  means  the 
binary  perceptron  did  not  converge  within  104  steps.  Even  though  a  convergence  proof 
does  not  exist,  preliminary  computer  simuations  do  indicate  a  bipolar  solution  weight 
pattern  can  be  obtained  provided  the  dimensionality  ,V  is  high  enough.  The  statistical 
capacity  with  N  =  16  is  shown  in  Figure  III. 3. 2.  Note  that  the  capacity  curve  starts  to 
roll  off  at  approximately  M,  the  number  of  stored  vectors,  equals  to  2/V/3,  whereas  the 
capacity  curve  of  a  linear  perceptron  (not  shown)  does  not  roll  off  until  \f  N .  This  can 
be  explained  by  the  fact  that  the  solution  region  of  a  linear  perceptron  may  not  contain  a 
bipolar  point.  In  such  case,  a  bipolar  vector  does  not  exist  which  classifies  all  input  vectors 
correctly. 


The  second  scheme  that  can  be  used  to  recognize  multiple  patterns  is  described  below. 
For  ease  of  comparison,  we  first  analyse  the  ordinary  (non-  bipolar)  linear  filter  of  the 
second  scheme  designed  for  the  desired  purpose.  Consider  the  following  algorithm:  the 
operation  to  be  performed  is  given  as 


r 


sgn 


‘.J 


j  +1  if  fm  belongs  to  vF; 
[  —  1  otherwise, 


(III. 3 .4) 


where  *F  is  the  class  of  objects  to  be  recognized.  The  filter  h(i,j)  is  generated  by  forming 
a  weighted  sum  of  all  the  input  patterns,  i.e. 


M 

^  {1 1 1.3.00.) 

m' 


where 


f  +1  if  fmi  belongs  to 'F; 
l—l  otherwise. 


(II  I.3.ob) 


To  see  that  the  above  scheme  is  capable  of  recognizing  multiple  objects,  rewrite  equation 
(III  3.4)  as 


sgn 


N-  M 


O, 


r  ^  ^  j)fm'  (*>  j) 


i ,]  m' 


{III.  3.6) 


The  RHS  of  the  above  equation  is  composed  of  two  terms,  namely,  the  signal  (first)  term 
and  the  noise  (second)  term.  Provided  that  the  signal  term  is  sufficiently  large  comparing 
to  the  noise  term,  the  correct  response  is  expected. 

Our  next  task  is  to  obtain  a  theoretical  estimate  of  the  capacity  of  the  filter.  YVe  define 
the  capacity  to  be  the  number  of  vectors  Xfc  that  can  be  stored  in  a  filter  with  vanishingly 
small  probability  of  error  for  sufficiently  large  dimensionality  .V.  Assuming  fm(i,j)  to  be 
a  discrete  sequence  of  bipolar-valued  independent  random  variable,  i.e. 


P  fm(iJ)  -  \  P  -  11  2; 

<  (1 1 1.3.7) 

Pfm{iJ)  fm'WJ')  P'fmiiJ)'. 


it  can  be  shown  that 


•V/ 


.V2 

8 ln\  ’ 


(III.  3.8) 


Shift  invariance  can  be  incoporated  into  this  system  by  modifying  the  operation  to  be 


.v-' 


g(i',j')  =-  sgn  Y^h(i,j)fm(i  -t-  i',j  f  j')  -  {  ‘  J 


if  fm  belongs  to  'F; 
otherwise. 


(///.:?.<)) 


In  such  case,  the  capacity  Mc  can  be  shown  to  be 


Mc  s 


/V2 


16  InN 


{III.  3.10) 


The  digitally  computed  correlations  of  random  sequences  (statistics  given  by  eq. (111.3,7)) 
and  the  linear  weighed-sum  filter  with  different  number  of  vectors  stored  are  shown  in 
Figure  III. 3. 3.  It  can  be  seen  from  the  simulated  results  that  as  the  number  of  vectors 
stored  in  a  filter  increases,  so  do  the  magnitudes  of  the  sidelobec.  This  phenomenon 
accounts  for  the  limit  of  the  number  of  vectors  that  can  be  stored  in  a  weighed  sum  filter. 
The  histogram  of  capacity  with  dimensionality  N  equals  128  is  also  shown  in  Figure  III. 3. 4. 
The  mean  of  the  histogram  agrees  with  the  theoretical  result. 

For  the  bipolar  analog  of  the  above  scheme,  consider  the  same  operation  given  by 
equation  (III. 3. 4).  In  this  case,  the  filter  is  generated  by  bipolarizing  the  filter  given  by 
equation  (III. 3. 5).  It  is  given  as  follows, 


r  M 


h{i,j)  =  sgn 


(///.  3.11) 


Using  the  same  statistics  for  the  stored  vectors,  it  can  be  shown  that  the  capacity  of  a 
bipolar  weighted  sum  filter  is  only  reduced  by  a  factor  of  7r/2.Thus 


Mb 


2  yv2 

7r  1 6lnN  ' 


for  filters  with  shift  invariance  incoporated,  and 


Mb 


2  ;V2 
7r  8 InN' 


{III.  3.12) 


(// /.3. 13) 


for  filters  without.  Computer  simulated  correlations  of  random  vectors  and  the  bipolar 
filter  formed  by  the  above  algorithm  are  shown  in  Figure  III. 3. 5.  Note  that  the  magnitude 
of  the  peak  decays  as  the  number  of  vectors  stored  increases  in  this  case.  This  phenomenon 
is  understandable,  since  as  M  increases,  the  relative  information  of  each  vector  stored  in 
the  filter  decreases.  The  histogram  of  capacity  of  the  corresponding  bipolar  filter  is  also 
shown  in  Figure  III. 3.6.  The  mean  of  the  histogram  agrees  with  the  theoretical  result. 


III. 4  Optical  System  Implementation 

Elsewhere  we  have  described  a  two-dimensional  incoherent  acousto-optic  image  correlator5 
and  have  discussed  various  advantages  associated  with  this  form  of  processing.  Here  we 
will  describe  the  implementation  of  bipolar  correlations  using  such  a  system.  The  basic 
system  architecture  is  shown  in  figure  (III. 4.1)  and  consists  of  a  16  element  vertical  LED 
array  in  the  reference  plane,  an  acousto  -  optic  device  (AOD)  in  the  input  plane,  and  a 
CCD  camera  in  the  output  plane.  An  IBM  PC  is  used  to  load  reference  image  data  into 
the  electronic  memory  shown.  The  data  is  read  out  to  the  LEDs,  32  bits  at  a  time  in  such 
a  way  as  to  modulate  the  intensity  of  the  ith  LED  with  consecutive  pixel  values  of  the  ith 
reference  image  line.  A  2  bit  D-A  converter  preceeding  each  LED  allows  for  representation 
of  up  to  4  distinct  intensity  levels  per  pixel.  In  the  applications  reported  here,  three  of 
these  levels  are  used  to  represent  an  unbiased  bipolar  signal  (-1,0,1)  as  a  biased  unipolar 
signal  (0,1,2).  The  imaging  optics  and  signal  timings  are  such  that  the  TDI  output  of  the 
CCD  is  proportional  to  the  correlation  between  the  bipolar  input  and  reference  images 
plus  some  bias  term.  The  ith  LED  is  intensity  modulated  by  the  current  : 

At(t)  =  [1  +  a,(t)]rect(t/T’i) 

where  a,(t)  represents  the  bipolar  signal  associated  with  the  ith  line  of  the  reference  image, 
and  T\  is  the  reference  image  width.  The  jth  line  of  the  input  image  results  in  an  intensity 
modulation  of  the  diffracted  light  after  the  AOD  (assuming  uniform  incident  intensity) 
that  is  given  by  : 

Bj{t  -  x/v)  =  l  +•  bj(t  -  x/v)\2rect(t—^) 
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where  b}(t)  is  the  bipolar  signal  associated  with  the  jth  line  of  the  input  image,  and  T-2  is 
the  TV  horizontal  line  time  plus  blanking  interval.  Using  this  we  obtain  the  nth  line  of 
the  CCD  output  as  : 


Cn(x) 


n'  —  n  rOC 

y  /  An’ ^M-n{t)Bn>{t  +  x/v)dt 

'  =  n-  .Vf-  l  ^ 
n'=n  . T 

V  /  |1  *  an*+..vf_ri(t))[l  +-  bn>(t  4-  x/v)]2dt 

n'  -  rc  —  ,Vf  —  l  ^  T 

n'  —  n 

\{x)  t  2  ^  /  an>  +  M  n[t)bn>(t  +-  x/v)dt 

n'  n  \  J  ~T 


(///.LI) 


Tn(x) 


+-  xjv)dt 


where  Tn(x )  is  the  bias  signal  on  the  nth  output  line,  M  is  the  number  of  CCD  lines  and 
the  correlation  time  T  is  given  by  :  T  --  T\  +  7’2. 

In  order  to  obtain  the  true  correlation  a(x,n)  *b(x,n)  we  clearly  must  remove  the  bias 
term  from  equation  (III. 4.1).  It  is  important  to  note  that  under  the  condition  : 


where  the  symbol  <>  represents  the  expectation  or  average  value,  the  bias  signal  7'rl(x)  is 
independent  of  the  input  and  reference  signals. 

111. 5  Correlation  Results 

We  have  implemented  the  above  system  and  have  obtained  correlation  results  as  shown 
in  figure  (III. 5.1).  In  our  experiments,  the  three  images  shown  in  figure  (111.5.1a)  were 
presented  to  the  TV  camera  as  the  input  scene.  The  input  scene  signal  was  first  passed 
through  a  DC  block  in  order  to  remove  any  bias  .  After  also  removing  the  horizontal  and 
vertical  sync  pulses  the  signal  was  amplified  and  a  constant  level  was  added.  This  level 
was  adjusted  so  as  to  insure  that  the  input  signal  was  strictly  positive.  For  a  bipolar  input 
scene  then  ,  this  unipolar  signal  is  simply  the  [1  +  6(<)j  described  earlier.  This  modified 
signal  was  used  to  amplitude  modulate  a  50  MHz  RF  carrier  which  was  then  applied  to 
the  AOD. 

Various  reference  images  were  generated  using  an  IBM  PC.  Each  reference  image  was 
16  lines  high  by  32  pixels  wide  however,  by  virtue  of  the  bipolar  encoding  scheme  employed, 
a  data  array  of  32  bits  by  32  bits  was  necessary  to  represent  one  image.  At  the  beginning 
of  each  input  image  horizontal  line  time,  the  32  lines  of  reference  data  were  read  out  to 
the  LEDs.  Preceeding  each  LED  was  a  2  bit  resistive  ladder  network  used  to  convert  two 
bits  of  digital  data  into  one  of  the  3  appropriate  analog  levels.  The  output  correlations 
are  shown  in  figures  (III. 5. lb,c,d)  for  various  reference  images.  As  can  be  seen  from  the 
figures,  the  system  exhibits  good  cross  correlation  suppression  and  strong  autocorrelation 
peaks. 

111. 6  Bias  Considerations 

In  order  to  represent  a  bipolar  signal  in  light  intensity,  a  constant  bias  is  added  to  the 
signal  before  modulation.  The  presence  of  this  bias  term  in  A(t)  and  B(t)  results  in  the 
bias  term  Tn(x)  in  the  system  output.  Removal  of  Tn(x)  is  therefore  required  before  the 
true  correlation  may  be  observed.  Shown  in  figure  (III. 6.1)  is  a  typical  line  of  the  bias 
signal  Tn(x).  As  expected,  the  bias  signal  is  approximately  triangular,  characteristic  of 
the  correlation  between  the  two  rectangular  functions  associated  with  A(t)  and  B{t). 

Removal  of  the  bias  function  Tn{x)  is  achieved  by  first  generating  a  bias  reference 
image  and  a  bias  input  image.  A  recording  of  the  2-D  output  correlation  of  these  two 
images  is  exactly  Tn(x)  assuming  the  conditions  of  eq.  III.  1.2  .  Thereafter,  in  order  to 
retrieve  the  true  correlation  from  the  system  output,  7'n(x)  may  be  subtracted  from  Cn(x) 
on  a  pixel  by  pixel  basis.  This  will  yield  the  desired  signal  : 

a(x.  n)  r  b(x,  n)  Cn{x)  Tn(x) 

If  the  condition  of  eq.  III. 4.2  is  not  satisfied  then  the  bias  term  T„  (x)  becomes  signal 
dependent.  In  general,  it  is  easily  shown  that  the  bias  signal  is  given  by  : 

n'.Z?  rT  t  v 

Tn{x)  ^  /  2  -  2 n(t)  t  brif(t  x,v)rect{t  T\)rect(  )dt 

n'  -  n  .Vf  t- 1  '  T  - 

(///.«.  1  ) 
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We  see  that  the  first  term  is  the  typical  triangular  bias  ;  however,  the  second  term  is  in 
fact  signal  dependent.  Various  schemes  for  removing  the  second  term  from  eq.  III. 6.1 
may  be  envisioned.  One  particularly  straightforward  method  would  be  to  simply  record 
the  two  signals  : 

n'  =  n  rT  t  _  / 

signal  1=  V]  /  an.  _n(t)rect{t /Tx)rect( 1  )dt 

n'  =  n-  M  ±1  J  1  2 


n'  =  n  fT  t  _  x/v 

signal2  =  I  bni(t  -  x/v)rect(t/Tl)rect( - 1 )dt 

n>=n-MHJ~T  72 


and  perform  a  pixel  by  pixel  subtraction  off  line  as  we  did  with  the  signal  independent 
term.  It  is  important  to  note  that  if  the  quantity  of  interest  is  the  inner  product  of  two 
images  rather  than  the  entire  2-D  correlation,  then  the  condition  of  eq.  [III. 4. 2]  may  be 
relaxed  somewhat.  In  this  case  it  is  necessary  only  to  have  equal  numbers  of  Is  and  -Is 
over  the  entire  image  to  insure  signal  independent  bias. 

III. 7  Implementation  of  a  Linear  Descriminant  Function 


The  above  system  has  been  utilized  in  a  pattern  recognition  scenario.  Taking  advantage 
of  the  flexibility  afforded  by  the  computer  generated  reference  images,  it  is  possible  to 
form  arbitrary  linear  combinations  of  these  images  to  generate  more  powerful  filters.  A 
LDF  based  on  three  images  was  generated  using  a  bipolar  perceptron  type  algorithm  as 
described  in  section  III. 3.  In  our  experiments  the  algorithm  was  init : alized  with  the  filter  : 

•»'  =  £  lMG‘~  H  IMG , 

i  in  't'l  i  in  ^2 


where  'F t  and  ^2  are  the  two  classes  and  I MGt  is  one  of  the  three  images  to  be  classified. 
Since  only  three  levels  may  be  represented  in  our  system,  a  threshold  arithmetic  must  be 
implemented.  That  is,  when  executing  the  above  algorithm  we  have  that  1-1  1  and 

that  -l  -  1  =  -1.  Although  there  is  no  guarentee  of  convergence  in  the  bipolar  case,  we 
observed  convergence  for  all  three  dichotomies.  The  LDF  as  found  above  was  used  as  the 
reference  image  while  the  three  seed  images  were  placed  in  the  input  plane.  The  results 
for  three  different  class  assignments  are  shown  in  figure  (III. 7.1).  We  see  that  successful 
classification  was  achieved  for  all  three  of  the  nontrivial  dichotomies  possible. 


III. 8  Conclusion 

We  have  shown  that  the  performance,  in  terms  of  output  SXR  and  storage  capacity,  of 
bipolar  correlation  filters  does  not  severly  degrade  when  compared  to  conventional  cor¬ 
relation  filters.  As  a  matter  of  fact,  bipolar  SLMs  are  free  of  the  problems  of  limited 
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dynamic  range  and  non-linearity.  We  have  also  presented  experimental  results  pertaining 
to  an  incoherent  acousto-optic  correlation  architecture.  The  particular  system  we  describe 
implements  bipolar  correlations  between  a  TV  scene  and  a  computer  generated  reference 
image.  Good  correlation  results  have  been  obtained.  We  have  also  successfully  imple¬ 
mented  a  two  class  decision  machine  using  this  correlator.  Many  issues  still  remain  to  be 
researched.  For  instance,  how  a  bipolar  filter  is  optimized  is  one  of  them.  How  a  bipolar 
perceptron  is  trained  to  recognize  multiple  objects  is  another  example.  For  a  certain  set  of 
data,  a  bipolar  filter  which  will  make  the  right  decision  in  recognition  may  not  exist  at  all. 
It  is,  nevertheless,  the  authors’  belief  that  the  probability  of  having  such  a  bipolar  filter 
approaches  one  for  sufficiently  high  dimensionality. 
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Figure  III. 1.1.  Output  vs  input  SNRs  for  the  binary  and  conventional  correlators. 
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Figure  III. 2.1.  (d)  Image  of  of  the  letters  A,BiC  and  D.  (e)  Computer  simulated  correlation 
result  of  (a)  and  (d). 
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>.l.  The  rate  of  convergence  of  a  bipolar  perceptron.  The  rate  of  convergence  of 
perceptron  of  the  same  data  set  is  also  plotted  for  comparison  .  The  circles  on 
2  the  bipolar  vectors  generated  by  taking  the  sgn[  |  of  the  ordinary  perceptron 
orni/.e  the  data  set  correctly. 
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Figure  III. 3. 3.  The  simulated  correlations  of  the  linear  weighed-sum  filters  of  different  num¬ 
ber  of  stored  vectors  and  some  of  their  stored  vectors.  The  threshold  level  for  recognition 
is  one  half  of  the  energy  of  each  vector. 
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Figure  III. 3. 5.  Computer  simulated  correlations  of  bipolar  weighed-sum  filters  of  differ¬ 
ent  number  of  stored  vectors  and  some  of  their  stored  vectors.  The  threshold  level  for 
recognition  is  one  half  of  the  energy  of  each  vector. 
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Figure  III. 6.1  :  Typical  Line  of  Correlator  Bias 


Figure  III. 7.1  :  Discriminant  Function  Implementation 


(a)  Correlator  Output  for  Class  I  =  IMG1  +  IMG2  /  Class  2 

(b)  Class  1  =  IMG1  +  H  /  Class  2  -=  IMG2 

(c)  Class  1  =  IMG2  +  H  /  Class  2  =  IMG  l 


IV  MULTIPLICATIVE  TIME  AND  SPACE  INTEGRATING 
ACOUSTO-OPTIC  SPECTRUM  ANALYZER 

IV. 0  Introduction 

In  this  section  we  will  describe  multiplicative  time  and  space  integrating  (TSI)  acousto¬ 
optic  architectures  for  spectrum  analysis  of  images  and  of  I-D  signals.  The  TSI  approach 
is  used  which  combines  the  best  features  of  1-D  space  integrating  and  1-D  time  integrating 
spectrum  analyzers1.  Wagner  and  Psaltis  have  previously  demonstrated  experimentally 
an  additive  architecture  for  folded  spectrum  processing  of  I-D  signals2.  In  this  section 
we  describe  and  experimentally  demonstrate  a  multiplicative  processor  suitable  for  finding 
the  2-D  Fourier  spectrum  of  images.  The  information  signal  is  fed  into  a  Bragg  cell  and  a 
lens  takes  the  Fourier  transform  in  the  coarse  frequency  (x)  direction  of  the  image.  In  the 
second  Bragg  cell,  the  DFT  signal  is  entered  which  along  with  the  time  integrating  CCD 
calculates  the  Fourier  transform  in  the  fine  frequency  (y)  direction  of  the  image.  We  begin 
the  section  by  describing  the  principles  of  TSI  processing.  The  space  integrating  spectrum 
analyzer  is  described  using  1-D  time  signals  and  2-D  video  images  as  inputs.  Then,  the 
Discrete  Fourier  transform  (DFT)  based  time  integrating  spectrum  analyzer  is  described 
and  experimental  results  are  given.  The  possible  TSI  architectures  are  given  and  analysis 
is  carried  out.  Experimental  data  is  presented  and  system  performance  issues  such  as  bias 
removal  techniques  are  discussed.  We  end  with  a  summary  of  the  processor  and  its  possible 
applications. 

IV. I  Principles  of  TSI  Processing 


Fig.IV.l.l  shows  the  basic  setup  for  a  TSI  architecture  for  spectrum  analysis.  Bragg  cells 
are  used  as  input  SLMs  because  of  their  high  bandwidth.  The  input  signal  s(t)  is  a  long 
1-D  signal  that  cannot  be  enclosed  in  the  time  window  of  the  Bragg  cell.  This  leads  to 
representing  s(t)  as  a  space-time  raster  f(x,n)  where  n  is  the  laser  diode  pulse  number. 
This  raster  then  passes  through  a  Fourier  transforming  lens  with  power  in  the  x-direction 
giving  another  space-time  raster  Fz(u.n)  which  consists  of  spatial  Fourier  transforms  of  the 
windowed  signal  s(t).  ,\ext,  the  raster  Fr(u,n)  is  spectrally  resolved  in  temporal  frequency 
by  using  time  integrating  Fourier  transform  modules.  The  output  of  the  system  is  a 
folded  spectrum  of  the  input  s(t).  The  motivation  for  dividing  the  2-1)  Fourier  transform 
operation  into  a  space  integrating  lens  Fourier  transform  and  N  time  integrating  Fourier 
transforms  comes  from  the  separability  of  the  2-D  Fourier  kernel  as  noted  in  equation 
(IV. 1.1).  Figure  IV. 1.2  shows  the  basic  principle  of  folded  spectrum  systems  originally 
shown  by  Thomas3.  Possible  inputs  s(t)  to  a  folded  spectrum  system  can  be  a  narrowband 
highly  coherent  signal  hurried  in  white  noise  (SETI)  or  a  video  signal  from  a  camera 
pointed  at  an  image  f(x,y).2  A  video  signal  has  a  natural  rastering  mechanism, therfore  tin' 

38 


2-D  folded  spectrum  output  is  the  2-1)  Fourier  spectrum  of  the  input  video  image. 


n  =  .V  -  1 


—  •  *  l  /• 

F(u,v)  =  /  f(x,n)e~ }2w,iZdx  e 


(/V  Ml] 


IV. 2  Space  Integrating  Spectrum  Analyzer 

Figure  IV. 2. L  shows  the  optical  setup  for  the  SI  Spectrum  Analyzer.  The  field  amplitude 
at  ^he  detector  plane  is 

F{x',t)  =  f  s{t  -  -)e-j2wz^dxe-^t  =  [vS{uv)e-}2wuvt  *  Asinc(Au)  e^^1 
J  A  v 

where  v  is  the  acoustic  signal  velocity  in  the  Bragg  cell  and  i(t)  is  the  complex  diffracted 
signal. 

Using  the  reference  signal  r  for  hetrodyne  detection,  the  intensity  on  the  detector 
plane  is 


For  a  single  tone  input 


the  intensity  is 


I(x',t)  =  l'vS(uv)e~}2,rvut  *  B(u)  +  r! 


s(t)  =  acos{2nf't) 


I(x\t)  =  a  2 B2(u  -  — )  —  jr  2  2  a  rB(u  -  ~)cos(2n f’ t  Q) 

v  v 

Taking  the  pulsing  action  of  the  laser  diode  into  account  every  T  seconds,  the  intensity  on 
the  nth  pulse  is 

[(x',n)  =  a  2 B2(u  -  — )  4-  r  2  +  2  a  rBiu  -  —)cos(2ir6nT  +  fl)  (/V'72.1) 

v  v 

where  B(u)  is  the  general  blurr  spot  taking  into  account  the  finite  aperture  of  the  Bragg 
cell  and  acoustic  apodization  effects.  Also, 

/'  =■  f,  -  6 


a  --  a  eJ 

Kquation  (Iv.2.1)  shows  that  for  a  single  tone  input,  the  output  of  the  processor  con¬ 
sists  of  a  constant  bias  term,  a  signal  dependent  bias  term  positioned  according  to  the 
signal  coarse  frequency  /,.  and  a  sampled  temporal  modulation  at  the  signal  fine  fre¬ 
quency.  The  third  term  contains  the  original  input  signal,  except  it  has  been  heterodyned 
to  baseband  by  the  pulsing  action  of  the  laser  diode.  It  is  this  temporal  modulation  that 
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allows  us  to  use  time  domain  spectral  analysis  schemes,  such  as  the  DFT  algorithm  to 
resolve  the  signal  into  higher  resolution  spectral  components.  In  order  to  avoid  aliasing 
effects  between  the  bandlimited  signals  in  each  coarse  frequency  blurr  spot. we  need  to 
satisfy  the  Nyquist  criteria  .i.e 

1  2 


where  T  is  the  sampling  time  of  the  laser  diode  and  Ta  is  the  aperture  of  the  Bragg  cell. 
When  s(t)  is  a  video  signal  where  /„  represents  the  nth  video  line,  the  field  amplitude 
on  the  detector  plane  can  be  written  as 


E(u,n)zzK  Fn(u)  *  B(u)\ 


{IV.  2.2) 


where  K  is  a  scaling  factor.  For  a  camera  looking  at  a  tilted  grating  image  with  a  spatial 
frequency  Uo  in  the  x  coordinate  direction  (fast  temporal  variation)  and  frequency  vq  in 
the  y  (slow)  direction, 

/n(x)  =  cos[2nuo{x  +  n<j>)  ] 


where  Tv  is  the  video  line  time.  In  this  mode  of  system  operation,  the  laser  diode  pulsing 
rate  should  be  set  equal  to  the  video  line  rate. i.e.  T  =  Tv  For  example,  if  we  have  a 
grating  with  variation  only  in  the  x  direction,  (  vq  =  0),  from  equation  (IV. 2. 2)  we  get 


E{u,n)  =  K  B[u  -  uo)  +  B(u  +  uo) 


(/e.2.3) 


and  for  a  variation  only  in  the  y  direction  of  the  image, 


E(u.  n)xzK B(u)cos(27Tv0nX) 


(/V.2,1) 


X  =  vT 


From  equation  (IV. 2. 3)  it  is  clear  that  the  fast  spatial  frequency  variations  along  the  x 
direction  of  an  image  can  be  resolved  in  the  Fourier  plane  by  using  an  acousto-optic  device 
as  an  input  transducer  and  a  Fourier  transforming  lens  to  spatially  channelize  these  fast 
frequencies  into  coarse  frequency  bins.  Equation  (IV. 2. 4)  tells  us  that  spatial  variations 
of  frequency  t'o  along  the  y  direction  in  the  image  are  converted  by  the  pulsing  space 
integrating  acousto-optic  spectrum  analyzer  to  a  sampled  temporal  variation  of  the  same 
frequency  in  a  coarse  frequency  bin.  This  temporal  variation  can  then  be  channelized 
into  its  fine  frequency  bins  along  the  y  direction  by  using  the  DFT  algorithm  based  spa¬ 
tially  distributed  time  integrating  modules  generated  by  using  imaging  optics  and  a  C(.'l) 
detector  array.  Figure  IV. 2. 4  shows  the  grating  video  images. 


IV. 3  Time  Integrating  Spectrum  Analyzer  Experiment 


We  have  used  the  DFT  Algorithm  approach  to  time  spectrum  processing4.  The  DFT 
algorithm  is  described  by 


n  =  At  -  1 


s(k)  = 


n  —  0 


W,v  =  e  1  "v 


where  S(k )  is  the  1-D  Fourier  spectrum  of  the  input  signal  s(n).This  operation  can  be 
broken  up  into  a  summation  over  n  of  the  products  of  the  input  signal  samples  with  its 
corressponding  DFT  matrix  columns.  As  the  DFT  matrix  is  symmetric,  we  require  only 
half  the  terms  in  a  DFT  column.  The  analog  nature  of  each  column  is  sinusoidal  with  a 
half  cycle  increase  with  each  increasing  DFT  column.  The  DFT  matrix  represented  as  a 
1-D  signal  is  a  stepped  frequency  chirp  signal  that  can  be  written  as 


t  nT 

d(t)  =  ^2  cos[nAu>(t  —  nT)\rect[ — — — ]  (7V.3.1) 

71  =  0 

We  used  a  personal  computer  to  generate  the  DFT  signal.  The  DFT  mask  and  the  half  cycle 
step  behavior  are  shown  in  Figure  IV. 3.1.  This  method  calculates  the  Cosine  transform 
i.e.  the  real  part  of  the  Complex  Fourier  Spectrum.  The  spectrum  analysis  is  performed 
along  the  length  of  the  DFT  column  which  acts  as  a  spatially  distributed  local  oscillator 
(DLO)  in  time.  The  frequency  of  the  oscillators  increases  from  the  top  (DC)  of  the  column 
to  the  bottom  (half  the  Nyquist  sampling  frequency).  If  one  of  the  input  signal  frequencies 
coincides  with  one  of  the  DLO  beads,  it  beats  with  it  to  build  up  an  interference  peak 
at  the  location  of  that  DLO  bead.  In  this  way,  the  frequency  components  of  the  input 
signal  are  resolved  into  the  fine  frequency  bins  along  the  DFT  column  (y)  direction.  Figure 
IV. 3. 2  shows  the  optical  system  set  up  in  the  laboratory.  The  laser  diode  is  pulsed  at  the 
video  line  rate  and  the  input  signal  s(t)  into  the  Bragg  cell  is 

s(t)  =  cos[(u;o  +  <5)t]  +  d(t)cosu0t 


where  u>q  is  the  center  frequency  of  the  Bragg  cell.  The  signal  s(<)  is  imaged  on  to  the 
CCD  via  the  +1  diffracted  order  light  from  the  Bragg  cell  with  appropriate  single  sideband 
filtering  done  in  the  Fourier  plane  of  cylinderical  lens  C 2.  The  light  intensity  integrated 
on  the  CCD  is 
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)  and  3  are  signal  dependent  spatial  frequency  and  phase  terms  and  C  is  constant  bias 
term.  Typical  value  for  is  as  1.86  cycles  per  mm  which  gives  only  a  few  cycles  on  the 
CCD  surface  giving  a  very  low  modulation.  Therfore,  we  see  peaks  located  at  the  sine 
type  function  maximum  position  y  =  corressponding  to  the  signal  fine  frequency  6.  In 
our  experiment.  6  is  varied  from  DC  to  7.8  KHz  which  corressponds  to  the  bandwidth  of  a 
video  image  along  its  vertical  (y)  or  slower  direction.  Figure  IV. 3. 3  shows  the  experimental 
results  without  bias  removal.  The  peak  moves  along  the  imaged  DFT  signal  ridge  with 
changing  fine  frequency. 

IV. 4  Architectures  for  TSI  Spectrum  Analysis 

Figures  IV. 4.1  and  IV'. 4. 2  show  two  possible  TSI  architectures  for  spectrum  analysis 
using  crossed  Bragg  cells  ,  imaging  optics,  a  pulsing  laser  diode  ,  CCD  detector  array 
and  supplementary  electronics.  The  general  signal  chain  for  both  systems  is  as  follows. 
The  DFT  signal  d(t)  added  to  a  reference  bias  a  is  mixed  up  to  the  center  frequency  jjo 
of  the  Bragg  cell  before  being  used  as  input  to  the  acousto-optic  device.  This  composite 
signal  si(t)  is  then  imaged  on  to  the  CCD  plane  with  appropriate  single  sideband  spatial 
filtering.  In  the  other  orthogonal  Bragg  cell, we  enter  the  sum  of  a  reference  chirp  c(t)  with 
the  information  signal  v(t)  (eg.  video, tones)  that  has  been  mixed  to  the  Bragg  cell  center 
frequency.  A  lens  takes  the  Fourier  transform  of  the  signal  s2(t )  in  the  AO  device  and  this 
spectrum  is  made  to  coincide  with  the  detector  plane.  The  electric  field  incident  on  the 
CCD  detector  plane  is 

E  =  f  sl(t-~)e-^^xdx\{s2(t  +  y~)\ 

Ja  V  V 

where 

Si{t)  =  y(t)  —  c(t)'cosur)t 
s2(t)  =  | a  +■  d(t)  cosjjot 

n  =  N  -  1 

v(0  =  Y,  ~  nT)rect(t~Y~) 

n  -  O 

for  a  video  signal. 

n  =•  V  1 

Cit)  Yl  COSb(t~nT)2:reCt(t^fL~) 

n  1 » 

and  d(t)  is  defined  by  IV. 3.1.  Taking  into  account  the  pulsing  of  the  laser  diode  (<5(<  m.7’)) 

and  the  single  sideband  filtering,  we  get  the  intensity  on  the  nth  pulse  as 

E„  2  %  v ( f ( u v )  •  C{vu))  *  AsincAu]  a  e  e  Jn  ^  y  \  (/V.l.l) 

On  expanding  equation  ( / V'.  4.1)  (see  appendix  A),  and  collecting  similar  terms,  we  get 
the  total  time  integrated  interferomertically  generated  charge  distribution  recorded  on  the 


CCD  after  a  frame  time  to  he  given  by 


n  —  .y  -  i  ^ 

I(u,  y)  —  K Re  ^  (F„(u)  -  B(u))cos(n  y)\  ~  Ki  -  K2  (IvA.2) 

n  -0 

where  the  first  term  gives  the  2-D  Cosine  transform  of  the  input  signal  and  K\  and  K2  rep¬ 
resent  constant  bias  and  signal  dependent  bias  terms  respectively.  K  is  a  scaling  constant 
and  we  have  assumed  that  the  reference  chirp  spectrum  is  uniform  over  the  bandwidth  of 
the  information  signal.  The  general  blurr  spot  B(u)  has  been  used  in  the  analysis.  For 
a  continous  time  signal  as  input,  the  first  term  of  equation  (7V.4.1)  represents  the  folded 
spectrum  of  the  input.  On  the  other  hand,  for  a  video  signal  input,  it  gives  the  2-D  video 
image  spectrum.  See  appendix  A  for  the  impulse  response  of  the  system  corressponding 
to  the  two  types  of  input.  The  architecture  in  figure  IV. 4.1  although  analytically  sound 
has  two  important  practical  limitations.  In  order  to  achieve  high  spectral  resolution  using 
present  day  hardware  .i.e.AOD’s,  CCD’s,  the  focal  length  of  the  imaging  lenses  runs  into 
meters.  Secondly,  for  ideal  operation  the  single  sideband  filtering  has  to  be  done  in  the 
plane  of  the  second  Bragg  cell  which  is  not  possible  because  of  the  thickness  of  the  crystal 
and  its  support  casing.  The  architecture  in  figure  IV. 4. 2  solves  both  these  problems  by 
using  crossed  cylinders  for  imaging.  This  imaging  method  gives  lens  values  that  assure 
overall  short  system  length  and  desired  demagnifications  for  compatibility  with  the  CCD 
detector  specifications.  In  our  case,  we  used  10  cm  and  15  cm  cylinderical  lenses  to  achieve 
the  desired  demagnifications. 

Figure  IV. 4. 3  shows  the  laboratory  setup  of  the  architecture  in  Figure  IV. 4. 2.  The 
system  timing  is  controlled  by  the  CCD  detector  7.16  MHz  internal  oscillator.  The  CCD 
blanking  signal  is  used  to  generate  a  15.734  KHz  laser  diode  trigger  signal  with  a  pulse 
width  of  100  nsec.  This  corressponds  to  a  light  flash  every  video  line  time.  The  anamorphic 
gaussian  beam  profile  of  the  laser  diode  is  collimated  by  the  spherical  lens  such  that  the 
long  axis  is  along  the  x  direction  of  the  SI  AOD  aperture.  The  Bragg  diffracted  light 
from  the  first  AOD  is  Fourier  transformed  by  a  spherical  lens  in  the  x  direction  and 
recollimated  in  the  y  direction  giving  a  slit  of  light  positioned  in  the  aperture  of  the  second 
AOD.  The  Bragg  diffracted  light  from  the  second  AOD  is  imaged  on  to  the  CCD  using  the 
crossed  cylinders.  In  the  Fourier  plane  of  the  cylinder  with  power  in  the  y  direction,  the 
appropriate  spatial  single  sideband  filtering  is  done.  The  AOD's  used  in  this  experiment  are 
slow  shear  mode  Teliriurn  Oxide  devices  with  an  aperture  time  of  70  nsec  and  bandwidth 
of  10  MHz.  A  60  MHz  center  frequency  is  used  for  the  AOD’s.  The  reference  chirp  is 
stored  in  a  digitally  programmable  read  only  memory  (PROM)  and  read  out  each  laser 
diode  trigger  pulse  using  the  CCD  pixel  clock  and  a  digital-to-analog  converter.  Figure 
IV. 4. 4  shows  a  reference  chirp  oscilloscope  trace  that  was  used  in  the  experiments.  The 
chirp,  being  digitally  generated  is  perfectly  coherent  at  every  laser  diode  pulse.  The  I") FT 
signal  is  generated  by  an  IB.M-PC  image  processing  work  station  that  is  locked  to  the  CCD 
detector  clock.  For  video  signal  inputs,  the  CCD  camera  is  locked  to  the  CCD  detector 
array  by  using  special  synchronizing  circuitry.  For  folded  spectrum  processing  of  tones, 
the  laser  diode  trigger  frequency  is  adjusted  according  to  the  Nyquist  criteria. 


Experimental  results  from  this  processor  are  shown  in  figures  IV. -1.5  and  IV.  1.6  for  a 
fine  frequency  analysis  bandwidth  of  7.86  KHz.  The  processor  is  operated  with  a  coarse 
frequency  resolution  of  80  KHz  per  pixel.  Data  in  figure  IV. 4. 5  corressponds  to  a  fine 
frequency  variation  from  DC  to  7.8  KHz  with  a  zero  coarse  frequency  variation.  We 
see  the  interference  peak  move  from  left(DC)  to  right (7.86  KHz)  along  the  DC  coarse 
frequency  bin.  Figure  IV. 4. 6(a)  shows  the  2-D  folded  spectrum  of  a  single  tone  input 
where  the  position  of  the  peak  gives  the  coarse  and  fine  frequencies.  Coarse  DC  position 
corressponds  to  the  base  of  the  picture.  Figure  IV. 4. 6(b)  shows  the  2-D  video  spectrum 
of  a  tilted  grating  used  as  the  video  image.  The  DC  of  the  fast  (coarse)  variation  along 
the  x  direction  of  the  image  corressponds  to  the  central  bin  between  the  two  peaks  on  the 
pk  *ire.  The  slow  (fine)  variation  along  the  y  direction  of  the  image  corressponds  to  the 
left  to  right  movement  of  the  peak  along  each  fast  frequency  bin.  We  get  two  peaks  instead 
of  one  because  the  video  signal  is  mixed  with  60  MHz  before  entering  the  Bragg  cell  giving 
an  amplitude  modulated  signal  whose  spectrum  has  two  sidebands  located  symmetrically 
around  the  60  MHz  carrier  frequency.  The  single  tone  experiments  are  carried  out  using  a 
signal  generator  whose  frequency  is  varied  around  60  MHz.  The  2-D  results  corresspond 
to  spectrums  after  electronic  bias  subtraction. 


IV . 5  Bias  Removal  Tec h tuques 


Bias  separation  from  the  desired  signal  is  always  an  important  issue  in  interferomet  ric 
time  integrating  processors.  We  discuss  three  ways  in  which  the  bias  can  be  removed. 
Pixel  by  pixel  electronic  bias  subtraction  can  be  done  by  using  two  synchonously  operated 
CCD’s,  one  to  record  the  signal  plus  bias  terms  and  the  other  for  recording  bias  terms. 
In  our  experiments,  we  used  a  PC  based  image  processing  work  station  with  a  frame 
grabber  to  store  the  bias  for  subtraction.  Another  approach  is  to  place  the  desired  signal 
on  a  spatial  carrier  which  is  converted  to  a  temporal  variation  by  the  natural  readout 
mechanism  of  the  CCD  and  then  bandpass  filtered  by  appropriate  electronics.  In  the  TSI 
processor,  the  2-D  spectrum  can  be  placed  on  a  spatial  carrier  by  mixing  the  DFT  signal 
by  an  frequency  offsetted  Bragg  cell  center  frequency  signal  eg.  60.1  MHz  and  using  a  time 
delayed  reference  chirp  signal.  The  bias  subtraction  and  carrier  demodulation  methods  of 
bias  removal  do  not  solve  the  problem  of  limited  system  dynamic  range  as  processing  is 
done  after  the  detection  of  the  signal.  The  effective  system  dynamic  range  at  the  output 


-  DR  (f  rSBft) 


where  DR  is  the  dynamic  range  of  the  output  detector  (CCD)  and  SBR  is  the  signal  to 
bias  ratio  on  the  detector.  In  most  TI  systems,  SBR  is  much  less  than  one.  One  way  of 
solving  this  dynamic  range  problem  is  using  a  photorefractive  crystal  for  bias  removal1. 
The  desired  spectrum  is  generated  on  a  spatial  carrier  along  the  crystal  side  and  is  read 
out  by  a  Bragg  matched  beam.  The  possible  architecture  for  the  TSI  spectrum  analyzer 
using  a  photorefractive  crystal  is  shown  in  figure  IV.5.1.  This  architecture  has  definite 
attractions  such  as  higher  light  efficiency,  better  system  dynamic  range  and  no  reference 
chirp  requirement  in  the  signal  chain.  In  addition,  the  single  sideband  filtering  does  not 
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have  to  he  done  in  the  Bragg  cell  aperture  and  this  system  can  process  higher  space- 
bandwidth  signals  as  the  resolution  of  the  photorefractive  crystal  is  much  greater  than 
that  of  the  CCD  detector.  As  this  system  is  not  common  path,  it  is  more  sensitive  to 
jitter. 

IV. 6  Processor  Performance  Issues 

The  Caltech  processor  shown  in  figure  IV. 4. 3  achieved  coarse  and  fine  frequency  res¬ 
olutions  of  80  KHz  and  100  Hz  respectively.  The  number  of  coarse  and  fine  frequency 
resolution  bins  corressponded  to  384  and  491  respectively.  The  bandwidth  of  the  system 
is  controlled  by  the  bandwidth  of  the  reference  chirp  which  is  programmed  in  a  high  speed 
PROM  to  5  MHz  for  video  signal  inputs.  Modulation  depths  of  30  percent  were  obtained. 
The  Caltech  processor  performance  is  limited  by  nonuniformities  in  the  reference  chirp 
spectrum  due  to  inhouse  electronics,  the  limited  temporal  coherence  of  the  high  power 
hitachi  laser  diode  and  intermodulation  product  terms  resulting  from  AOD  and  amplifier 
nonlinearities.  The  coarse  frequency  resolution  is  limited  by  CCD  pixel  size,  focal  length 
of  Fourier  transforming  lens,  AO  cell  aperture  and  apodization  effects.  For  the  video  pro¬ 
cessor,  the  ideal  coarse  frequency  resolution  of  15.7  KHz  could  be  obtained  by  using  a 
larger  focal  length  lens  and  smaller  CCD  pixels.  The  fine  frequency  resolution  is  limited 
by  integration  time  of  the  CCD  which  corressponds  to  60  Hz.  Note  that  using  carrier 
demodulation  for  bias  removal  reduces  the  number  of  frequency  bins  along  the  carrier  di¬ 
rection.  For  using  a  photorefractive  crystal  such  as  Bismuth  Silicon  Oxide  (BSO)  for  bias 
removal  in  our  proposed  architecture  we  need  around  8  VV  pulsed  or  12.6  mW  continous 
power  to  record  gratings  with  an  Argon  laser  beam  3.  This  high  power  requirement  can 
be  reduced  for  test  purposes  if  we  wait  a  few  video  frames  before  reading  the  integrated 
spectrum  on  the  crystal  using  the  Helium-Neon  readout  beam.  The  principle  of  operation 
involves  the  coherent  addition  of  the  image  spectrum  over  successive  video  frames  to  build 
up  enough  charge  density  for  recording  gratings  on  the  crystal.  Another  important  issue  to 
note  is  that  in  a  multi-tone  signal  input  environment,  the  TS1  processor  will  not  introduce 
crosstalk  because  each  tone  is  spatially  separated  in  the  Fourier  plane  thus  giving  cross 
terms  that  have  negligible  contributions. 


Appendix 

This  appendix  gives  the  derivation  of  equation  [V.4.2  starting  with  equation  IV. 4.1 
given  earlier. For  simplicity  of  the  analysis,  we  have  assumed  the  u  -  0  axis  to  corresspond 
to  the  position  of  the  Bragg  cell  center  frequency  (u,’0)  on  the  CCD  plane.  The  intensity 
on  the  nth  pulse  given  in  equation  IV. 4. 1  can  be  written  as 

In  =  iS  ~C)(a±d)2  (1) 

where  '5'  and  'C'  are  Fourier  spectrums  of  the  information  signal  and  reference  chirps 
respectively  taken  by  the  space  integrating  Bragg  cell  spectrum  analyzer  and  ‘a’  and  ‘<f’ 
are  the  reference  bias  signal  and  the  spatially  filtered  and  hetrodyned  to  baseband  DFT 
signal  respectively,  ‘a’  is  a  real  quantity  and  ‘d’  is  given  by 

d  =  e~]n^Ly 

d\2  =  1 

Using  the  above  relations  and  expanding  equation  (1)  we  get 

In  =  2Re[aSC'Re{d)\  +  (1  +  a2),C\ 2  +  2\CfRe{ad) 

-  (L  -  a2)|S|2  +  2\S\2Re{ad)  +2(1  +  a2)Re(SC ' )  ^ 

Note  that  the  reference  chirp  is  a  perfectly  coherent  chirp  with  no  dependence  on  n 
as  the  same  digital  chirp  is  generated  every  laser  diode  pulse.  Re  in  equation  (2)  stands 
for  real  value  of  a  function.  Equation  (2)  consists  of  6  terms  where  the  first  term  gives 
the  2-D  spectrum,  the  second  and  third  terms  give  the  constant  bias  term  and  finally  the 
fourth, fifth  and  sixth  terms  give  the  signal  dependent  bias  term.  Assuming  the  reference 
chirp  spectrum  C  is  uniform  over  the  information  signal  bandwidth  and  has  an  arbitary 
fixed  phase  which  can  be  taken  to  be  zero  for  simplicity  of  analysis,  we  can  write  C  - Cn 
where  Co  is  a  constant.  Now  taking  the  summation  over  all  the  .V  laser  diode  pulses  to 
find  the  charge  distribution  on  the  CCD  ,  the  first  term  in  equation  (2)  gives 

n  =  N-l 

h=2aC0Re[  ]T  S Re(d)'  (3) 

n  —0 

Substituting  in  equation  (3)  the  value  of  S  for  a  video  signal  and  the  value  of  d  .  we  get 
the  desired  2-D  spectrum 


n—S— 1  , 

/,(u.y)  KRe  (/'„(«)  *  0(u))cos( H  y)| 

*——*  v 
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where  K  is  a  constant  term  involving  the  constant  terms  Cq  and  a  and  the  finite  size  of 
the  CCD  detector  array  in  the  u  and  y  directions.  We  have  used  the  generalized  blurr  spot. 
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£?(u)  to  account  for  the  finite  aperture  of  the  Bragg  cell  in  the  SI  direction.  The  second 
and  third  terms  of  equation  (2)  on  summation  give  the  constant  bias  term  h\  to  be 

„  ^  2,,  o  m  2  ,  .sin(i\3y)  ... 


K i  =  A'Co  (1  a2)  -  2clCq  cos(ay) — (5) 

sin(3y) 
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where  the  first  term  in  equation  (5)  gives  a  constant  bias  level  and  the  second  term  gives  a 
fine  frequency  DC  ridge  along  y  —  0  with  a  very  slow  spatial  modulation  undetectable  by 
the  CCD.  Now  taking  the  summation  over  n  of  the  fourth, fifth  and  sixth  terms  in  equation 
(2),  we  get  the  signal  dependent  bias  term  K 2  as 


K2  =  iV(l  -r  a2)|S(u)|2  +  2a|S(u)j2cos(ay) 


sin(Npy) 


2C0(1  +  a2)|5(u)a1 


"■4  —  -  / 1  ^  v  —  / 1  — sin[i3y )  u'  1 

Note  that  the  absolute  value  operation  on  S(u,n)  in  equation  (2)  removes  the  time  depen¬ 
dence  in  the  spatial  Fourier  transform  taken  by  the  lens  giving  only  the  coarse  frequency 
bias  ridge  terms  in  equation  (6).  The  first  and  third  terms  in  equation  (6)  give  bias  terms 
positioned  at  the  signal  coarse  frequency  ridges.  The  constant  a\  depends  on  the  fine 
frequency  of  the  signal.  The  second  term  gives  the  coarse  frequency  ridge  crossed  with  the 
fine  frequency  DC  ridge.  Combining  the  results  from  equations  (4),  (5)  and  (6)  we  get  the 
total  time  integrated  charge  distribution  on  the  CCD  to  be 

n  —  jV  -  l 

[(u,y)^KRe  ^  {Fn(u)  «  B(u))cos(n-~y))  +  Ki  +  K2  (7) 

n  < ) 

The  impulse  response  of  the  system  for  a  video  image  input  with  a  spatial  frequency  un  in 
the  x-direction  and  spatial  frequency  v,,  in  the  y-direction  of  the  image  is 

I[^CO.s(r])  ’J  B(u  Un)  -  B(u  r  Un)]  (8) 

n  (.V  1) 

.  y 


where  we  used  the  fact  that  the  modulation  along  y  is  slow.  We  will  get  two  peaks 
positioned  at  the  coordinates  u  Un.y  “7r,'n^  and  u  -  Un,  y  -  27rt'u^^. 

For  a  single  tone  input  signal  with  /,.  and  ff  the  coarse  and  fine  frequencies  respec¬ 
tively.  we  get  the  impulse  response  of  the  system  as 

li(U-  )  (9) 

sin[  r 


nffT  - 


This  time  we  get  one  peak  located  at  u 
dependent  constant  phase  term. 


*ffT  aL,  vvhcrc  s  is  a  fine  frequency 
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V.  Adaptive  Broadband  Array  Processing 


V.O  Introduction 

Adaptive  signal  processing  techniques  for  processing  broadband  phased  arrays  are  pre¬ 
sented.  After  discussing  the  essential  differences  between  narrowband  arrays  and  broad¬ 
band  arrays,  optimal  forms  for  the  broadband  space-time  filters  are  derived.  Adaptive 
processing  systems  which  converge  to  the  optimal  forms  are  then  presented.  Acousto-optic 
implementations  using  multi-channel  acousto-optic  devices  are  presented.  One  particu¬ 
lar  implementation  makes  use  of  a  photorefractive  crystal  as  a  time-integrating  detector 
[11,121. 


V.l  Narrowband  Processors  in  Broadband  Environments 

The  extension  of  adaptive  filtering  techniques  to  the  space-time  domain  is  applicable 
for  sonar  and  radar  signal  processing  where  the  outputs  from  an  array  of  sensors  must  be 
weighted  and  summed  to  optimally  estimate  a  signal  in  the  presence  of  noise  '  1-3] .  When 
the  array  processor  encounters  only  narrowband  signals  centered  at  a  common  carrier 
frequency,  /0,  the  filter  structure  shown  in  Fig.  V.l  is  identical  to  that  of  the  temporal 
filter  except  that  the  input  samples  are  the  outputs  of  the  sensors  instead  of  a  tapped 
delay  line.  In  such  applications,  only  spatial  discrimination  of  signals  is  required,  and  the 
jammer  nulling  capabilities  of  such  systems  are  quite  good  as  long  as  the  jammers  are 
narrowband.  The  output  can  be  described  by  the  following  equation 

v 

40  -  (VM) 

n  -  l 

where  .V  is  the  number  of  elements  in  the  array  and  hn  is  the  adaptively  controlled  weight 
for  the  nth  sensor  output  un(t). 

After  suppressing  the  temporal  carrier  term  exp{j2n fct},  the  signal  received  by  the 
nth  element  due  to  a  narrowband  signal  arriving  at  an  angle  8  from  boresight  can  be 
expressed  as 

un(t)  -  a(t)exp{  j'2n f^.sinOnd/c},  ( V '.‘2 ) 

where  c  is  the  speed  of  light,  a (t)  is  the  slowly  varying  envelope  of  the  signal,  and  d  is  the 
spacing  between  the  array  elements.  The  goal  of  the  adaptation  is  the  minimization  of 
the  mean  squared  error  between  the  actual  output  of  the  processor  and  a  desired  signal. 
The  scenario  with  which  we  will  be  concerned  is  the  one  where  the  desired  signal  is  known 
to  be  on  boresight,  and  interference  from  directions  other  than  this  must  be  suppressed. 
This  is  sometimes  known  as  sidelobe  cancellation  1  .  An  example  of  narrowband  nulling 
is  shown  in  Fig.  V.2,  where  the  array  gain  pattern  for  a  uniformly  weighted  array  and  the 


nm  iuiii! 


minimum  mean-squared  error  gain  pattern  for  receiving  a  signal  on  horesight  and  rejecting 
a  jammer  of  the  same  frequency  at  a  sidelobe  angle  are  plotted. 

Although  the  jammer  is  received  at  a  sidelobe  angle  and  is  attenuated  relative  to  the 
signal,  the  power  of  the  jamming  interferences  is  typically  much  larger  than  that  of  the 
signal,  and  so  the  summed  output  of  the  uniformly  weighted  processor  will  be  dominated 
by  the  interference.  In  this  case,  the  optimum  pattern  to  which  the  adaptive  processor 
converges  places  a  perfect  null  precisely  at  the  direction  of  the  interference. 

Given  the  size  of  the  array  (.V  -  the  number  of  receiving  elements),  an  upper  bound 
on  the  number  ot  spatial  jammers  of  the  same  frequency  /.,  that  the  array  is  able  to  reject 
can  be  established.  The  mth  jammer  from  a  set  of  Xf  jammers  can  be  represented  by  an 
.Y-element  vector  u("l)  whose  nth  element  is  given  by 

=exp{-j2nfllsinOmndc}.  (l'.3) 


Simultaneous  nulling  of  all  of  the  XI  jammers  leads  to  the  system  of  equations 

/  u(!1}  t41}  ...  u|v1}  \ 


(M) 
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V  *. 


,(M) 
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\  kN  / 

(V.4) 


where  s  represents  the  signal  vector,  and  the  last  equation  of  the  system  arises  from  the 
constraint  that  the  signal  be  accepted.  If  the  jammers  are  distributed  spatially  such  that 
the  vectors  and  s  are  linearly  independent,  then  in  order  for  a  solution  h 

to  exist,  there  can  be  at  most  .Y  equations  in  the  system,  yielding  the  upper  bound  of 
\f  <  .X  1.  The  linear  independence  condition  is  akin  to  that  of  general  position  in 
pattern  recognition  9  and  assures  that  the  jammers  are  irregularly  distributed  in  space. 

Since  the  narrowband  arrays  assign  only  one  weight  per  sensor,  no  consideration  of  the 
temporal  content  of  the  received  signal  is  taken.  The  spatial  distinction  of  a  signal  is  judged 
based  on  the  observed  differences  in  the  phase  of  the  signal  as  received  by  each  element. 
If  a  particular  jammer  contains  many  frequencies  spread  over  a  significant  bandwidth, 
each  component  will  result  in  a  unique  relative  phase  difference  from  element  to  element 
and  hence  will  look  like  a  multiplicity  of  jammers  directed  at  different  angles.  Thus,  if  ,i 
narrowband  processor  encounters  a  sufficiently  wideband  interference  in  some  direction, 
it  will  use  up  all  of  its  degrees  of  freedom  (the  adjustable  weights)  to  null  this  single 
directional  noise. 

Two  examples  are  shown  in  Figs.  Y.3  and  V.l.  where  the  minimum  mean-squared 
error  array  patterns  are  plotted  for  a  boresight  desired  signal  and  a  multi-frequency  jammer 
incident  at.  0  60  .  The  element  spacing  in  each  case  is  ,,'y  .  The  jamming  frequencies  are 

.85/.,.  and  .7 with  the  signal  power  and  the  power  in  each  jamming  frequency  being 
equal.  For  each  case  (Fig.  V.3,  Fig.  V.l),  a  gain  pattern  for  each  jamming  frequency  is 
plotted,  since  the  array  response  is  frequency  dependent.  Fig.  V.3  depicts  the  response  of 
a  2  element  array,  and  it  is  clear  that  only  the  frequency  .85 /0  is  nulled  perfectly,  while  Fig. 
V.l.  which  corresponds  to  a  1  element  array,  shows  perfect  nulling  for  all  three  frequencies. 
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GAIN  PATTERN  FOR  UNIFORM  WEIGHTS 
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ANGLE  FROM  BORESIGHT  (IN  DEGREES) 


Fig.  V.2  N  arrowband  Sidelobe  Cancellation  (arrow  indicates  the  an¬ 
gular  position  of  the  jammer):  a)uniformly  weighted  array  gain  pat¬ 
tern;  b)minirimm  mean  squared  error  array  gain  pattern  (jammer  inci 
dence  angle  is  50 


2  ELEMENTS-  f=.85f0 


-60  -30 


-60  -30 


ANGLE  FROM  BORESIGHT  (IN  DEGREES) 


ANGLE  FROM  BORESIGHT  (IN  DEGREES) 


2  ELEMENTS  f = . 7f 0 


-60  -30  0  30  60  9' 

ANGLE  FROM  BORESIGHT  (IN  DEGREES) 


Fig.  V.3  Minimum  Mean  Squared  Error  (Iain  Patterns  for  a  Multiple 
Frequency  Jammer  (/  -  /0,  .85/,,,  .7 /,,)  Incident  at  CO  from  Bore- 
sight  (2  element  array,  element  spacing  is  c/2/,,):  a)response  for  /,,, 
b)response  for  ,85/o,  c)responsc  for  .7 f0 


4  ELEMENTS  f-fO 


1.2 


4  ELEMENTS  f«.85f0 


Fig.  V.4  Minimum  Mean  Squared  F>ror  Gain  Patterns  for  a  Multiple 
Frequency  Jammer  (/  /q,  .85/n,  . 7 /, , )  Incident  at  60  from  Bore- 

sight  (4  element  array,  element  spacing  is  c/2/0)'.  a)response  for  /o, 
b)rcsponse  for  .85/o,  c)response  for  .7 f0 
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V.2  Optimum  Broadband  Systems 


In  narrowband  applications,  bandpass  filters  are  used  to  assure  a  narrowband  signal 
and  noise  environment  as  shown  in  Fig.  V\l,  but  where  signals  with  large  bandwidths 
must  be  processed,  such  filters  cannot  be  used.  Even  in  broadband  noise  environments, 
we  would  like  the  X'  element  array  to  be  capable  of  cancelling  .V  -  1  jammers  in  general, 
regardless  of  their  respective  bandwidths.  This  requires  that  more  degrees  of  freedom 
be  available  for  the  output  of  each  sensor  than  the  single  weight  that  the  narrowband 
processor  provides.  Shown  in  Fig.  V.5  is  an  .V  element  array  processor  that  satisfies 
the  requirements  for  operation  in  broadband  environments  by  passing  the  output  of  each 
sensor  through  a  linear  time-invariant  filter,  or  equivalently,  a  tapped  delay  line  5-7  . 
While  spatial  discrimination  is  still  made  possible  by  the  spatial  sampling  done  by  the 
array  of  sensors,  the  additional  capability  of  making  distinctions  based  on  the  temporal 
content  of  the  received  signals  is  offered  by  the  array  of  filters. 

An  optimum  choice  of  the  linear  filters  based  on  the  mean  squared  error  criterion  can 
be  made  in  the  present  scenario  of  a  desired  signal  incident  on  boresight  and  broadband 
jammers  from  other  directions.  Let  s[t)  represent  the  desired  signal  which  is  assumed 
to  derive  from  a  stationary  random  process  with  zero  mean  and  autocorrelation  function 
given  by  R(t)  =  E{s(t  +  r)s'(t)j.  Since  the  desired  portion  of  the  received  signal  arrives 
on  boresight,  there  is  no  dispersion,  and  each  element  receives  an  identical  desired  signal 
as(f).  The  interference  that  arrives  from  different  directions  will  induce  signals  in  each 
element  with  different  delays. 

Let  i'n(t)  be  the  interference  component  received  by  the  nth  element,  assumed  also 
to  derive  from  a  zero  mean  stationary  random  process  with  covariance  function  given  by 
~lmn(r)  =  E  vm(t  -  r)tr(0  .  The  total  signal  received  by  the  nth  element  is  therefore 
given  by  un(t)  —  o,s(f)  -  vn(t).  The  output  of  the  array  processor  is  given  by 


N  i-  oo 

(Q  -  y  ^  /  hn{r)un(t  -  7 -)dr. 

n  -  1  J  ~  00 


The  problem  is  to  minimize  the  mean-squared  error 

c  Mr)  -  E  z(t)  ■  s(t)  2’, 


(V.5) 


(l -.6) 


by  varying  the  filter  kernels  hn(r).  By  using  simple  variational  arguments,  it  is  shown 
in  Appendix  B  that  the  minimizing  solution  for  hn{r )  must  satisfy  the  system  of  linear 
integral  equations 


v  r  *> 

rri  1  ' 


hm(r)  nR[t  r)  ■  -)rnn(t  r);  dr  R(t). 


W. 7' 


where  n  1.2 . V.  'IVansformat iori  of  the  above  system  to  flip  frequency  domain  vie 

,Y 


Hm(f)  «5(/)  *  rm„(/)  S(f). 


(V.Ha) 


m  l 
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Fig.  V.5  Frost’s  Broadband  Array  Processor 
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Here  S(f)  is  the  spectral  density  function  of  the  signal  and  rm„(/)  is  the  cross-spectral 
density  matrix  of  the  interference,  given  by  the  Fourier  transforms  of  the  corresponding 
covariance  functions,  and  Hm{f)  is  the  frequency  response  of  the  rnth  filter  of  the  array. 
F'or  comparison,  the  corresponding  equation  for  a  narrowband  processor  with  a  single 
complex  weight  wrn  for  the  rnth  receiving  element  is 

v 

U'm  ot  -  Crnn  .  =  1,  for  n  =  1.2 . N,  (V.Sb) 

tn  -  1 

where  a  is  the  amplitude  of  the  received  signal  component  on  boresight  and  Cmn  - 
E  vmv'n  ,  the  covariance  matrix  of  the  noise  components.  A  comparison  of  Eqs.  Y.8a 
and  V.8b  shows  that  the  optimal  broadband  strategy  is  to  simply  provide  an  optimum 
narrowband  weight  for  each  frequency  /.  Thus,  the  optimum  broadband  system  reduces 
to  a  channelized  system  of  optimum  narrowband  sub-systems. 

If  we  now  assume  that  the  signal  spectrum  is  approximately  white  (e.g.,  spread  spec¬ 
trum  codes)  with  spectral  density  So  and  the  received  signal  component  is  small  compared 
to  the  noise  so  that  aSo  <C  1,  then  Eq.  V.8  can  be  approximated  by 

N 

Y  Hm{f)rmn(f)  *  S„.  (Y.9) 

m=  l 

The  frequency-dependent  response  of  the  individual  filters  can  now  be  varied  to  compensate 
for  the  spectral  characteristics  of  the  environment,  whereas  the  narrowband  processor 
offered  no  such  flexibility  with  its  single,  frequency-independent  weight  for  each  sensor. 


V.3  Directional  Cancellation  of  a  Single  Broadband  Interference[2,4-6 


The  example  of  a  cancellation  of  a  multi-frequency  interference  given  in  Fig.  V.l 
shows  that  a  narrowband  array  processor  suppresses  a  broadband  jammer  by  placing  the 
null  over  a  wide  region  near  the  interference  direction.  This  is  true  in  general  for  any  nar¬ 
rowband  array  processor  with  a  number  of  elements  sufficiently  large  to  null  a  broadband 
interference.  A  large  number  is  needed  because  the  array  pattern  shifts  with  frequency. 
By  solving  Eq.  V.9  explicitly  for  the  case  of  a  single  broadband  interference  incident  at.  an 
angle  0,  with  respect  to  boresight.  we  will  show  analytically  that  the  broadband  system 
suppresses  the  interference  with  a  null  that  does  not  shift  with  frequency,  placed  precisely 
at  9,. 

Let  a(t)  be  the  broadband  interference  waveform  incident  at  an  angle  0t  with  respect 
to  boresight.  so  that  the  nlh  element  receives  the  relatively  delayed  version  n(t  riA).  where 
A  ilsinO.  r.  Assume  that  the  interference  waveform  is  derived  from  a  zero  mean  station¬ 
ary  random  process  with  the  autocorrelation  function  given  by  y(r)  E  a(t  •  r)n'(f)  . 
the  covariance  between  the  interference  received  by  the  nth  and  rnth  elements  is  given  by 


1  (r) 


E  a (/  •  r  mA)o‘(f  tiA) 
-,(r  -  (n  rn) A). 


( V 10) 
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The  inclusion  of  sensor  or  receiver  noise  that  is  uncorrelated  for  each  element  and  assumed 
to  be  white  with  spectral  density  o*  yields 


lmn(r)  =  Tf(r  -r  (n  -  m)A)  ~  Smn6(r)o'; 


(VM1) 


as  the  covariance  of  the  total  noise  present  in  the  rnth  and  nth  elements.  The  noise  cross 
spectral  density  matrix,  required  for  the  optimal  solution,  is  the  Fourier  transform  of  the 
above  expression.  It  is  given  by 


r mn (/)  =  exp{j2wf(n  -  m)A}Sa(/)  -  6mno‘ 


( V'.  1 2 ) 


The  particularly  simple  form  assumed  by  the  noise  cross  spectral  density  matrix  en¬ 
ables  the  inversion  of  Eq.  V.9  to  be  carried  out  explicitly.  It  can  be  shown  that  the  solution 
of  Eq.  Vr.9,  which  gives  the  optimum  form  for  each  filter,  is  given  by 


Hn{f)  -  (So/cr^) {  1  -  exp{  -j7r/A(iV  +  I  -  2 n)} 


sin(nfNA)„  , 


(V13) 


The  gain  pattern  that  results  from  this  choice  of  filters  can  be  determined  by  computing  the 
array  response  to  a  monochromatic  plane  wave  signal  at  various  frequencies  and  incidence 
angles.  Specifically,  if  we  let  9  represent  the  incidence  angle  of  this  probe  beam,  the  gain 
pattern  as  a  function  of  the  probe  frequency  and  angle  is  given  by 

r  a\  —  ( _ r  e<  \r  i  \  \  '\  c  /„ 2\  I  sin(nfN A  ) 


g(f,0)  =  (exp{  -  jnf(: V  -  l)A'}So/<^){  — 

t  sin 


(a/A') 


-  exp{  j  tt  f  (/\r  -  A)(iV  4-  0  -H) 

sin(Trf  A) 

■sinjnf N ( A '  -  A ) ) _ Sy (f)  i 

sin{xf{ A'  -  A))  (of-  NSv(f)j  J  ’ 

where  A'  dsind  c.  At  incidence  angles  other  than  that  of  the  interference.  0,,  the  second 
term  is  small  since  A  i-  A'  so  that 


g(f,B  *  9t)  (S0jo*)exp{  -jnf(N  d-  1)A'} 


n  sin( 7r /.V A') 


( V.  15) 


j  ~  i  ;  ~  i)  / - r  k  j  j  \  ‘  /  —  )  /  r  \  i\  \  *  ' ) 

Sin(7T  J  Jv) 

which  is  simply  the  array  pattern  with  uniform  weighting. 

When  we  probe  near  the  interference,  however,  so  that  9  «  9 ,,  then  the  second  term 
becomes  appreciably  large  and  the  gain  becomes 


g(f,9  ss0t)  =  (S.,  o,*)exp{  jxf(\  -  L) A} 


sin(nfS  A) 


'  J  '  '  sin  ( 7T  /A)  (V.  Ifi) 

{1  \S4f):(o*  -  .V.V, ,(/))}. 

With  the  assumption  that  the  interference  power  is  large  compared  to  that  of  the  detector 
noise  so  that  .V.S ’,  (/)  •  a*,  the  gain  becomes  zero  near  the  interference  direction.  .Note 

also  that  this  null  in  the  array  pattern  remains  fixed  for  all  frequencies  where  there  is 


V.4  Adaptive  Array  Processor 


ts 


As  in  the  strictly  temporal  case,  the  adaptive  broadband  sidelobe  canceller  also  uti¬ 
lizes  convolvers  and  correlators  to  accomplish  the  desired  task.  However,  because  of  the 
multi-dimensional  nature  of  the  signals  (spatial  and  temporal),  arrays  of  correlators  and 
convolvers  must  be  employed.  The  basic  system  is  shown  in  Fig.  V.6. 

The  output  of  the  processor  is  simply  the  sum  of  all  of  the  signals  from  the  array 
elements  and  a  feedback  signal,  which  is  derived  from  a  cascade  of  multi-channel  correlation 
and  convolution  operations  performed  on  the  output  and  input  signals.  The  output  is  fed 
back  to  the  first  block,  which  calculates  the  correlation  between  the  output  signal  and  the 
,V  input  signals  from  the  array  elements.  Each  of  the  resulting  correlation  functions,  given 
by  ( un(t )  *  z(t)),  n  =  1,2, ....  N ,  is  then  convolved  with  its  corresponding  input.  The  .V 
convolved  results  are  summed  to  give 


v 

U  =  G  U„(Q  *  (un(Q  »  g(0)  (V.17) 

n—  i 

as  the  feedback  signal,  where  G  is  the  feedback  gain. 

To  show  that  the  adaptive  processor  approximates  the  optimum  response,  however, 
the  output  must  first  be  expressed  in  the  form  of  Eq.  V.5,  where  the  impulse  respon  e 
for  each  array  filter  is  explicitly  shown.  We  will  show  that  the  equation  describing  the 

filter  impulse  response  functions,  hn(t),  n  -  1,2 . S,  of  the  adaptive  processor  is 

approximately  equivalent  to  that  of  the  minimum  mean-squared  error  processor  given  by 
Eq.  V.9. 

In  order  to  clarify  the  analysis  to  follow,  wre  redraw  the  schematic  of  Fig.  V.6a  as 
shown  in  Fig.  V.6b  so  that  as  in  the  Frost  architecture  (Fig.  V.5).  the  output  from 
individual  filters  can  be  identified  arid  described.  Here,  the  output  of  the  nth  filter  is  given 
by 

'Jn(t)  J  hn{t  r)un(r)dr . 

or  equivalently.  }'„(/)  ff  nif)l'n(f)  where  Yn(f)  is  the  Fourier  transform  of  i /„(/).  The 
overall  output  is  then  given  by  c(/)  V^r|  y,x(t).  Inspection  of  Fig.  V '.Ob  yields 
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Taking  the  Fourier  transform  of  the  above  equation  and  rearranging  yields 
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on 


X>nm  -GU'n{f)Um{f)\Hm{f)  =  I, 


(V.20) 


where  lrn{f)  is  the  Fourier  Transform  of  u„(<).  If  the  feedback  gain  G  is  sufficiently  large 
so  as  to  leave  the  first  term  of  Eq.  V.19  negligible  compared  to  the  second,  then  Eq.  V.20 
becomes 

(/)»„(/)  =  i-  (V-.21) 


Thv  is  of  the  same  form  as  the  optimum  equation,  Eq.  V.9.  provided  that  the  input  SNR 
is  sufficiently  low  so  that  we  can  identify  u‘n(f)um(f)  as  the  estimate  of  the  noise  cross 
spectral  density  matrix  rmn(/)-  The  particularly  simple  form  of  Eq.  V.20  can  be  inverted, 
however,  to  give  a  closed  form  solution  for  the  adapted  filtering  functions.  The  derivation 
is  analogous  to  that  leading  up  to  Eq.  V.13  and  yields 


*.(/)  =  .-  .  n  =  1,2, IV. 

i  +  c£L,HU/)lJ 


(V.22) 


Suppose  that  the  total  signal  received  by  the  array  consists  of  a  weak  probe  signal  on 
boresight  represented  by  p(t)  whose  Fourier  transform  is  P(f )  and  strong  interference 
signals  incident  in  other  directions,  the  noise  received  by  the  nth  element  represented 
by  vn(t)  whose  Fourier  transform  is  Vn(/).  Thus,  the  total  input  is  given  by  un[t)  -- 
p{t)  +  vn(t).  The  probe  is  sufficiently  weak  as  to  only  negligibly  affect  the  determination 
of  the  filtering  functions  Hn(f )  so  that  approximately, 


Hn[f)  *  1 
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n  =  1,2,.., /V. 


(V.23) 


The  output  is  described  by 


Z(f)  =  £>»(/KM/) 


K(f)  ELlVmif) 


(V.24) 


P(f)  *  Vn{f)\  • 


Further  expansion  yields 
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As  expected,  the  second  term  is  effectively  suppressed  if  the  feedback  gain  can  be  made 

large. 


V.5  Optical  Implementation 

The  optical  implementation  of  the  adaptive  sidelobe  canceller  is  a  fairly  straight¬ 
forward  extension  of  the  optical  techniques  used  for  the  temporal  systems  described  in 
Chapter  3.  The  extension  is  made  simply  by  using  arrays  of  convolvers  and  correlators 
to  handle  the  N  input  signals  from  the  antenna  array  sensors.  The  input  spatial  light 
modulators  to  be  used  are  multi-channel  and  single  channel  acousto-optic  devices  that 
offer  large  dynamic  range  and  can  operate  on  broadband  signals.  We  concentrate  first 
on  a  space-time  integrating  system  that  involves  the  use  of  a  photorefractive  crystal  to 
perform  the  time  integration.  A  strictly  space  integrating  array  processor  is  described  in 
Section  V.6. 

The  basic  task  of  the  processor  is  to  form  the  feedback  signal  given  by  Eq.  V.17. 
Note,  however,  that  since  correlation  can  be  expressed  as  un(t)  *  z(t)  —  un(  t)  *  ^(t),  and 
convolution  is  both  an  associative  and  commutative  operation,  the  feedback  signal  can  be 
expressed  equivalently  as 


,v 

“(0  -  G~(0  *  Y  un(0  * «;(  0-  (v'-2fi) 

n  -  1 

This  rearrangement  suggests  that  we  can  modify  the  system  architecture  somewhat  by  first 
computing  an  autocorrelation  function  of  each  input  signal,  adding  the  autocorrelations 
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and  convolving  the  summed  result  with  the  output,  ^(t),  to  form  the  feedback  signal  as 
illustrated  in  Fig.  V.7. 

For  the  purpose  of  discussion,  however,  we  adhere  to  the  strict  interpretation  of  Eq. 
V.17  as  discussed  in  V.4.  As  shown  in  Fig.  V.6a,  the  output  is  first  correlated  against 
each  of  the  X  input  waveforms,  and  each  resulting  correlation  is  then  convolved  with  its 
corresponding  input.  The  X  convolution-correlation  signals  thus  produced  are  summed 
to  form  the  feedback  signal.  The  multi-channel  correlations  are  performed  using  time 
integration  and  the  convolution  with  the  output  signal  is  achieved  with  space  integration. 

Shown  in  Fig.  V.8  is  the  multi-channel  correlator,  where  the  output  autocorrelations 
are  written  as  modulations  of  an  index  grating  formed  on  the  crystal  by  the  photorefractive 
effect.  More  specifically,  AODl  and  AOD2  are  arranged  to  operate  in  a  coherent,  additive 
architecture.  The  single  channel  AODl  diffracts  a  portion  of  its  input  light.  This  passes 
through  AOD2  without  being  affected  because  of  Bragg  mismatch  and  is  imaged  onto  the 
crystal  to  yield  the  amplitude  z(t  i-  x/v  -  T /2)exp{j2irfo(t  +  x/v)}  in  the  crystal  plane. 
The  multi-channel  AOD2  diffracts  the  undiffracted  beam  from  AODl,  which  is  aligned  to 
be  well  matched  to  its  acoustic  gratings  to  yield  the  optical  amplitude 

^  un[t  -  xjv  -  T/2)rn(y)exp{j2irfo{t  -  x/v)}, 


also  imaged  onto  the  crystal.  Here  rn(y )  characterizes  the  vertical  confinement  of  the 
acoustic  beam  in  each  channel  of  the  multi-channel  AOD  and  can  be  approximated  by 


rn{y)  i^fy)rect 


by 


n  =  1,2,..,  X, 


(V.27) 


6y  being  the  acoustic  beam  width  and  yo  the  separation  between  neighboring  channels, 
where  X  is  assumed  to  be  odd. 

These  two  amplitudes  are  incident  at  an  angle  with  respect  to  each  other.  Since 
both  have  the  same  Doppler-shift  in  frequency,  an  intensity  grating  is  seen  by  the  crystal. 
In  fact,  an  array  of  one-dimensional  modulated  index  gratings,  each  confined  vertically 
by  r„(y),  is  formed  within  the  crystal.  It  can  be  shown  that  the  modulation  functions 
are  proportional  to  the  desired  correlation  functions,  and  specifically,  the  diffracted  light 
due  to  the  nth  grating  when  the  crystal  is  illuminated  by  a  plane  wave  is  given  by  8 
En{x,y)exp{  -  j2nf,)(2x >)}.  Here 
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In  this  equation,  the  exponentially  decaying  window  function  has  been  approximated  l>\  a 
rectangular  window  of  temporal  duration  t.  The  resulting  integral  is  proportional  to  the 
correlation  function  of  un(t)  and  z(t)  with  2x  v  as  the  shift  variable. 


To  complete  the  task  of  forming  the  feedback  signal,  we  combine  the  previously  de¬ 
scribed  correlators  with  a  space  integrating  convolver,  resulting  in  the  system  shown  in 
Fig.  V.9.  The  added  portions  are  the  multi-channel  AOD  (AOD3)  which  is  driven  by 
the  system  output  z(t),  some  imaging  optics,  and  a  single  photodetector.  Since  the  auto¬ 
correlations  are  available  at  a  compressed  horizontal  spatial  scale,  2 x/v,  a  combination  of 
cylindrical  and  spherical  lenses  Cl,L5,L6  is  used  to  anamorphically  image  AOD3  onto  the 
crystal.  This  anamorphic  imaging  provides  a  2:1  demagnification  ratio  in  the  horizontal 
direction  and  1:1  in  the  vertical. 

A  portion  of  the  He-Ne  read  beam  is  diffracted  by  AOD3  and  passes  through  the 
crystal  unaffected  since  its  propagation  direction  is  not  properly  Bragg-matched  to  the 
correlation  gratings  in  the  crystal.  The  DC  beam  from  AOD3,  however,  is  arranged  to 
be  Bragg  matched  to  the  gratings  and  reads  out  the  correlation  functions.  If  the  acoustic 
velocity  in  AOD3  is  v,  the  same  as  that  of  AODl  and  AOD2,  and  all  of  the  AODs  are 
driven  at  the  same  frequency,  /o,  then  the  read  beam  diffracted  by  the  crystal  and  that 
diffracted  by  AOD3  are  collinear  and  interfere  temporally  at  /o  because  of  the  Doppler  shift 
induced  by  AOD3.  The  diffracted  light  amplitude  just  behind  the  crystal  due  to  the  crystal 
gratings  is  proportional  to  En(x,  y)exp{j2nf0(-2x/v)},  and  the  diffracted  light  from 
AOD3  evaluated  at  the  same  plane  is  exp[j2nf0(t  -  2x/v)]un(t  -  2x/v  -  T/2)rn(j/). 

The  sum  of  the  above  two  amplitudes  are  Fourier  transformed  by  lens  L7,  and  the 
resulting  intensity  distribution  is  integrated  across  the  Fourier  plane  by  a  sufficiently  large 
photodetector.  If  w  is  the  width  of  each  AOD,  the  result  is  the  output  photocurrent 
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un(t  -  2x/v  —  T /2)dt'dx }  +  DCterms. 

The  high  frequency  term  centered  at  fo  is  approximately  equal  to  the  desired  feedback 
signal  given  by  Eq.  V.17.  The  output  of  the  system  z(t )  is  then  formed  by  subtracting  the 
feedback  signal  from  the  sum  of  the  input  signals  Vn  un(t)  to  yield 
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un(t  2d  -  T/2)dt'd0\, 


where  J  x  v,  T  -  W/V,  and  G  is  the  feedback  gain.  Here  we  have  assumed  that  the 
signals  received  by  the  array  are  all  centered  at  f0  and  that  they  drive  the  AODs  directly 
after  being  amplified. 


By  defining  an  equivalent  impulse  response,  hn(a ),  for  the  system  such  that  z(t) 
V ]n  J //°  hn(a)un(t  -  a)da ,  Eq.  V.30  can  be  expanded  further  to  give 
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(V.31) 


An  equation  that  describes  the  impulse  response  can  then  be  derived  by  noting  that  the 
output  appearing  on  the  left  side  of  the  above  equation  can  be  expressed,  using  the  impulse 
response.  With  the  assumption  that  the  integration  time,  r,  is  long  enough  to  warrant  the 
approximation  sinc[fr j  ss  6(f),  the  resulting  impulse  response  equation  is  given  by 


Mtf)  53  H&)  ~  (f)rect 
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hn(a)U;(f')Um(f')df'da, 


where  Un(f)  is  the  Fourier  transform  of  the  nth  input  signal  un(t).  Taking  the  Fourier 
transform  of  the  above  equation  gives 

»Af)  =  l-c£  r  -  f)T)sinc(f  -  f)T\Hm(f')VM')U„(ndf.  (V.33) 

m 

This  is  very  similar  to  the  optimal  equation  discussed  in  V.2  if  we  identify  the  product 
m  (/)  as  an  approximation  of  the  required  cross  spectral  density  matrix.  In  par¬ 
ticular,  the  effect  of  the  finite  convolution  time  is  to  distort  the  spectral  properties  of 
the  input  noise  field.  Thus,  the  optical  implementation  is  expected  to  place  a  spectrally 
broader  null  in  comparison  with  the  interference  bandwidth,  but  the  spatial  characteristics 
remain  similar  to  the  optimal  case. 


V.6  Adaptive  Array  Processor  with  Variable  Look  Direction 

The  array  processor  described  thus  far  is  a  sidelobe-cancelling  system  where  the  max¬ 
imum  sensitivity  (the  ‘‘look  direction”)  is  always  constrained  to  be  in  the  direction  of 
boresight.’  In  this  section,  an  acousto-optic  processor  is  presented  which  is  an  extension 
of  the  temporal  active  processor  described  in  to  the  space-time  domain.  The  optical  archi¬ 
tecture  considered  uses  space  integration  as  opposed  to  the  photorefractive  implementation 
considered  in  the  previous  section.  The  optimizing  criterion  used  is  that  of  maximizing 
the  output  S.\R  (signal  detection)  and  the  signal  need  not  arrive  on  boresight. 

The  scenario  considered  is  that  where  a  signal  waveform  .s(t)  is  incident  on  the  arra\  at 
a  known  angle  9  with  respect  to  boresight.  and  noise  (possibly  broadband)  from  directional 

*  The  look  direction  can  be  changed  by  introducing  appropriate  delays  in  the  signals 
received  by  each  element  before  they  are  processed. 


sources  corrupt  the  received  signal  waveform.  Let  sn(t)  =  s(t  -  n{dic)cosd)  be  the  signal 
waveform  received  by  the  nth  element.  The  total  signal  received  by  the  nth  element  is  then 
given  by  un(t)  =  sn(t)  -F  vn(t),  where  vn(t)  is  the  sum  of  all  noise  terms.  Here  d  is  the 
array  element  spacing,  and  c  is  the  speed  of  light. 

The  output  of  a  general  space-time  filter  can  be  expressed  as 


A  /-oo 

y(t)  =)j  un(T)hn(t  -  r)dT, 
n  =  l 


( V.3-1) 


where  un(t)  is  the  complex  envelope  of  the  total  rf  signal  (centered  at  fo)  received  by  the 
nth.  array  element,  and  hn(t)  is  the  filtering  function  for  each  channel.  A  similar  expression 
is  obtained  for  the  optically  implemented  space-time  filter  using  two  multi-channel  AODs 
shown  in  Fig.  V.10.  This  system  is  coherent,  and  the  output  can  be  shown  to  be  given  by 


/.\  2c ! 

y(0  =  ~y 


N  rT/2 

'-Zf 

n  =  lJ° 


un(t  -  t  -  T/2)hn(t  +  r)dr, 


(V.35) 


where  the  filtering  function  hn(t)  is  used  to  drive  the  second  AOD,  and  Cj  is  a  constant 
that  depends  on  laser  power,  AOD  diffraction  efficiencies,  and  the  quantum  efficiencies 
of  the  detectors.  The  only  differences  from  the  general  filter  (Eq.  V.34)  are  limited 
accumulation  time  and  the  time  compression  of  the  output.  This  is  of  little  consequence 
since  signal  detection  rather  than  estimation  is  considered  here.  The  noise  vn(t)  present 
in  the  received  signal  is  modeled  by  a  zero-mean  random  process  with  a  covariance  matrix 
given  by 

7mn(0  =  E[vm{t)v'n[t  -  r)j.  (V.36) 

It  can  be  shown  through  variational  arguments  that  the  choice  of  hn(t),  for  which  the  SNR 
of  the  AO  space-time  filter  is  maximum  at  a  specified  time  t0,  must  satisfy  the  following 
system  of  integral  equations: 
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The  details  of  this  derivation  are  given  in  Appendix  B. 

For  adaptivity,  we  need  to  calculate  and  continuously  update  the  filter  function  to 
drive  the  output  to  the  optimum  result.  As  with  the  temporal  active  processor,  the  output 
must  be  correlated  with  the  input  to  produce  the  appropriate  filter  function.  Since  the 
array  processor  has  .V  inputs  and  one  output,  this  requires  that  we  correlate  .V  signals 
with  a  common  one.  This  can  be  achieved  with  the  arrangement  shown  in  Fig.  V.l  1.  That 
figure  shows  the  use  of  a  multi-channel  AOD  driven  by  the  \  antenna  element  outputs  in 
conjunction  with  a  single  channel  AOD  driven  by  the  processor  output  signal. 

The  complex  envelope  of  the  output  of  the  nth  detector  element  is  given  by 


Fig.  V.10  M  ulti-Channel  Acousto-Optic  Space  Integrating  Convolver 


Fig.  V.ll  Multi-Channel  Acousto-Optic  Space  Integrating  Correlator 
with  a  Single  Reference 


where  c2  >s  a  constant  that  depends  on  the  laser  power.  AOl)  diffraction  efficiencies,  and 
the  quantum  efficiency  of  the  detector.  For  proper  correlation  to  appear  at  each  output, 
the  signal  driving  the  single  channel  AOD  must  be  time-compressed  by  a  factor  of  two. 
This  is  indeed  the  case  for  the  system  described.  Thus,  the  AOD  implemented  space-time 
filter  and  the  .V  channel  correlator  with  a  single  reference  are  compatible. 

Shown  in  Fig.  V.1‘2  is  the  array  processor  block  diagram  that  shows  the  interconnec¬ 
tions  required;  it  is  a  direct  extension  of  the  temporal  •>  'ive  processor  to  two  dimensions. 
The  output  from  each  antenna  element  is  correlated  w  ith  the  processor  output  to  produce 
the  filter  function  for  that  element.  The  “steering  vector"  $’(-()  determines  the  look 
direction  of  the  array  and  also  is  the  temporal  reference  signal  used  for  the  detection  of 
the  desired  signal  s (/). 

Fig.  V.13  shows  the  optically  implemented  adaptive  array  processor  with  the  AOD 
implemented  space-time  filter  in  the  upper  branch  and  the  .V  channel  correlator  in  the 
lower  one.  By  combining  Eqs.  V.35  and  \  .38,  the  equation  that  determines  the  filter 


function  hn(t)  is  seen  to  be 


( fi.39) 


where  (7  is  the  feedback  gain,  finder  conditions  of  low  input  SNR  and  large  feedback  gain. 
Eq.  A. 39  can  be  transformed  to  the  frequency  domain  to  yield 
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where  >'„(/)  is  the  Fourier  transform  of  .s„(t). 

For  comparison .  consider  the  Fourier  transform  of  Eq.  \  .36;  the  optimum  filter  equa¬ 
tion  is  given  b\ 

. (/)//-,(/)  A  erp(  j2n  f  t,,)Sn  (/),  (V.U) 

"i 

where  I  ,.  ,  (/(  is  the  spectral  « i » * ri ^ i t >  matrix.  Identifying  the  integral  in  Eq.  V.IO  as  the 
smoothed  estimate  of  the  cross  spectral  density  matrix  of  the  input  noise  vector,  Eq.  V .10 
is  approximately  equivalent  to  Eq.  V.ll.  However,  the  effect  of  the  finite  time  integration 
window  i-  -eon  m  the  -moot hing  ol  the  noise  spectrum. 
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Appendix 

Optimum  Broadband  Array  Processors 


The  optimum  broadband  array-processing  systems  discussed  in  Chapter  V  are  derived 
in  this  appendix.  The  two  optimality  criteria  considered  are  the  mean  squared  error  and 
the  maximum  output  SNR  conditions.  Simple  variational  arguments  are  used  to  derive 
the  necessity  conditions  for  these  criteria.  Since  the  costs  to  be  minimized  are  quadratic 
functionals  of  the  impulse  responses  that  are  varied,  these  necessity  conditions  are  sufficient 
as  well,  and  so  describe  the  systems  uniquely.  Throughout  this  appendix,  the  noise  received 
by  the  nth  array  element  is  represented  by  un(t)  which  is  assumed  to  be  a  zero  mean, 
stationary,  random  process  with  covariance  -y mn(r)  =  E[vm  (*K(*  -  01- 


Minimum  Mean  Squared  Error  Sidelobe  Canceller 

We  model  the  desired  signal  by  a  zero- mean  stationary  random  process  s(t)  whose 
autocorrelation  function  is  given  by  i?(r)  =  E[s(t)s'(t  -  r)].  Let  the  total  signal  (desired 
signal  plus  noise)  received  by  the  nth  element  be  represented  by  un(t)  =  s(t)  +vn(t).  Since 
the  scenario  of  interest  focuses  on  a  signal  arriving  on  boresight,  each  element  receives  the 
same  desired  signal  without  dispersion. 

The  form  of  the  processor  is  shown  in  Fig.  V.5,  where  the  parameters  that  must  be 
optimized  are  the  .V  impulse  response  functions.  The  output  signal  is  given  by 

N  r  OO 

z (0  ~Y2  l  hn{r)un(t  -  r)dr,  (.41) 

n=l  J~x 


and  the  cost  to  be  minimized  is  given  by 
C  hn(t)  E  z(t )  s(t)  2\ 

^2^2  f  [  hm(3)h»(T)  R(d  r)  *  lnm{3  r)  drdJ 

*  R(o)  iYl  J  Mr)#(  T)dr  •  r 


.42) 


where  we  have  assumed  the  noise  and  signal  portions  to  be  statistically  independent. 
Let  hr,(t)  represent  the  optimum  filter  for  the  n,h  element,  and  let  represent  any 

other  impulse  response  function.  In  order  For  h„[t)  to  be  the  cost-minimizing  solution, 
the  perturbed  cost  func  tion  ('  h„(t)  ■  mihl  have  a  minimum  value  at  the  point 

o  n h  ■  jit 1  0,  where  o h  and  0/  are  real.  The  two  conditions  that  express  tlm 

mathematically  arc1 

2  C  hjt)  nqjt )  0. 

’  1  (CM 

'  r  hjt)  ■  n./.Jt) 

fht  / 


0 


Carrying  out  the  above  calculations  leads  to  two  equations  which  combine  to  yield 


d3  -  R'(r)}dr  —  0.  (A4) 

Since  gn(t)  is  an  arbitrary  function  in  Eq.  A4,  the  terms  within  the  braces  must  sum 
identically  to  zero.  This  gives  the  following  condition  that  the  optimum  filter  must  satisfy: 

h-m(0)[R(/3  -r)+  lnm((3  -  r)\d(3  -  R‘(r)  =  0,  (A5) 

which  is  both  necessary  and  sufficient  because  the  cost  functional  is  quadratic.  This  is  the 
generalization  of  Wiener’s  result  to  ihe  space-time  domain  of  broadband  phased  arrays. 
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Maximum  Output  SNR  Array  Processor 

We  now  consider  the  task  of  signal  detection,  using  the  array  processor  of  Fig.  V.5. 
The  signal  is  a  known  waveform  s(t),  which  arrives  at  an  angle  9  with  respect  to  boresight, 
and  the  noise  field  is  the  same  as  that  discussed  for  the  MMSE  (Minimum  Mean-Squared- 
Error)  processor.  The  total  signal  received  by  the  nth  element  is  therefore  given  by 

“n(0  =  «n(0  ~  wn(0»  (A6) 

where  s„(t)  -  s(t  -  ndcosO/c)  is  the  signal  waveform  as  seen  by  the  nth  element. 

The  output  is  given  by 


e(f)  - 


r)hn(r)dr. 


(.47) 


Since  the  processor  is  linear,  the  signal  component  can  be  readily  identified  from  the  noise 
at  the  output..  The  output  signal  and  noise  are  described,  respectively  by 


^  -■  1  j  n  1 1  ( t) 

en,  , ..  (/) 


r)hn(r)dr. 


r)h„(r  )</,- 
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The  goal  is  to  maximize  the  output  SNR  given  by 


s  \  H 
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at  a  prescribed  t irne  t 


An  equivalent  problem  is  to  minimize  the  output  noise  power  while  constraining  all 
admissible  impulse  responses  to  give  the  same  output  signal  amplitude  at  t0.  In  particular, 
the  functional  expressed  by  Eq.  A9  can  be  maximized  by  minimizing  the  cost  given  by 


9SM0I  =  -Elkno.dc^o)!2!  -  /  5n^°  ~~  T)hn(T)dT  -  c} 

=  EE  /  f  h*{T)h'm{P)lnm{P  ~  T)drdp  (,410) 

n  m  J  J-°c 

-  /  Sn(to  -  T)hn(T)dT  -  c}. 
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where  the  constraint  estgnai(to)  =  c  has  been  included  through  the  use  of  the  Lagrangian 
multiplier  A. 

Again,  representing  the  optimum  filter  by  hn(t)  and  an  arbitrary  function  by  gn(t), 
we  apply  the  conditions  given  by  Eq.  A3  to  obtain 


r  oo  _  r  oo 

y!  /  /  "1 <nm{0  ~  T)hm(l3)d(3  -  Asn(t0  -  r))dr  =  0.  MU) 
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Since  gn{t)  is  an  arbitrary  function,  the  terms  within  the  braces  must  sum  identically  to 
zero,  yielding  the  condition 

E  /  ~lmn(r  .1)hm(j)d3  A  '  t).  (.412) 


which  is  the  matched  filter  result  generalized  for  phased-array  processing.  The  Lagrangian 
multiplier  A  is  immaterial  since  it  merely  multiplies  all  .V  filters  by  a  constant  factor  and 
has  no  effect  on  the  output  SNR.  Therefore,  Eq.  A 12  describes  the  optimum  result  to 
within  a  constant  factor 


